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Prólogo 


La computadora de bajo coste Timex/Sinclair 1000 suele adquirirse 
como una primera introducción al campo de las computadoras personales. 
Pero su aspecto de juguete defrauda a muchos presuntos compradores que 
no se percatan de su potencial considerable. Prueban unos pocos progra- 
mas de demostración, encuentran su uso dificultoso y le dejan a un lado. 
Este libro, escrito específicamente para el Timex/Sinclair 1000, abre nue- 
vos campos de aplicaciones para el estudiante universitario o de segunda 
enseñanza e incluso para el ingeniero y para el científico. 

Un estudiante universitario, o de segunda enseñanza, reconocerá mu- 
chos problemas aquí tratados, tales como la aritmética de los números 
complejos, combinatoria, distribuciones y soluciones a ecuaciones diferen- 
ciales y transcendentes, como los temas de sus asignaciones para trabajo en 
caso. Y los profesionales encontrarán que con estos programas pueden re- 
solverse fácilmente problemas matemáticos y de ingeniería con esta peque- 
ña computadora. 

Este libro contiene una colección de programas que personalmente he 
desarrollado durante varios años para su uso en mi profesión como inge- 
niero de comunicaciones y para mi entretenimiento en tareas de cálculos 
caseros y diseño electrónico. Abarcan una amplia gama de problemas en 
los campos de la ingeniería eléctrica, probabilidad, estadística, teoría de las 
colas, fiabilidad, ajuste de curvas, generación de gráficos, teoría de los nú- 
meros, informática, inteligencia artificial y otras materias. 

Algunos de estos programas se desarrollaron originalmente para las cal- 
culadoras programables HP-25 y HP-67 de Hewlett-Packard y se traduje- 
ron después al BASIC. En el proceso de traducción, añadí varias caracterís- 
ticas relativas a una mejor emisión de llamadas orientativas y de conten- 
ción de errores de entrada. Muchos más programas útiles se escribieron 
originalmente en BASIC y se añadieron a la colección. Unos pocos progra- 
mas se han publicado en revistas. En estos casos, se cita la revista de 
procedencia con los reconocimientos correspondientes. 

La codificación se ha desarrollado a partir de las versiones anteriores, 
menos sofisticadas, concebidas para las computadoras personales Radio 
Shack/Sharp. La traducción al BASIC del Timex/Sinclair 1000 exigió mu- 
chos cambios, pues los lenguajes BASIC de las computadoras Sharp y 
Timex/Sinclair sólo son parcialmente compatibles. Muchos programas se 
tuvieron que volver a escribir completamente y hubo de recurrirse al em- 
pleo de características funcionales avanzadas del Timex/Sinclair 1000 tales 
como la manipulación de cadenas y la visualización en pantalla. En el cur- 
so de la escritura de los programas del Timex/Sinclair 1000 desarrolló 
muchos artificios prácticos de programación que se detallan en el apéndice 
de este libro. Estos artificios, y diversas técnicas de programación descritas 


en la sección de Observaciones sobre la Programación de cada capítulo, de- 
ben aumentar las aptitudes y conocimientos prácticos de programación del 
lector. 

Todos los programas de este libro están especificamente escritos para la 
computadora con al menos, 2K de RAM o una computadora ZX-81 análo- 
gamente configurada. Sin embargo, cada listado va precedido por una 
descripción del programa y las ecuaciones fundamentales correspondien- 
tes, que deben facilitar la traducción a otras computadoras de BASIC. 
Cuando se traduzcan a otros lenguajes de computadoras, ha de tenerse cui- 
dado con las «idiosinerasias» de cada lenguaje. Por ejemplo, la función de 
enteros INT(- 3.5) proporcionará — 4.0 en BASIC, pero dará - 3.0 en For- 
tran (!). 

Espero que encuentre estos programas de utilidad y que le permitirán 
disfrutar todavía más de su Timex/Sinclair 1000. 


Cass Lewart, Holmdel, NJ, 1981 


Introducción 


Cada capítulo de programa, consta de cinco secciones: 

1. Descripción del programa. Proporciona al lector los conocimientos 
básicos necesarios y la comprensión del problema. 

2. Instrucciones. Ponen de manifiesto las secuencias claves requeridas 
para ejecutar el programa en la computadora Timex/Sinclair 1000. 

3. Ejemplos. Sirven para aclarar más el programa. Esta sección incluye 
párrafos sobre el enunciado del problema, discusión de los resultados y un 
ejemplo de listado del programa ejecutado directamente por la computado- 
ra por medio de la orden COPY. 

4. Observaciones sobre la programación. Explican las técnicas de pro- 
gramación y flujo del programa. También sustituyen a las sentencias de 
comentarios (REM) en los listados que se omiten debido a las limitaciones 
de memoria de la computadora de 2K. 

5. Listados de programas. Se imprimen con la orden LLIST y con in- 
terconexión de una impresora de matriz de puntos normal, correspondien- 
do cada listado exactamente al programa que genera los ejemplos. 

Observará que todos los programas-incluidos en este libro tienen algu- 
nas características comunes. Son fáciles de utilizar y proporcionan unas 
llamadas, o mensajes, de orientación y una contención de errores de carác- 
ter completo. Los programas no son triviales; v.g. no hacen algo que podría 
hacerse más fácilmente con una calculadora manual. Con muy pocas ex- 
cepciones, cada programa realiza, sin esfuerzo, funciones que serían difíci- 
les de efectuar de otra forma. Ninguno de los programas es simplemente 
una traducción de fórmulas. Y no utilizan trucos «indocumentados» para 
ahorrar un octeto (byte) de memoria a costa de claridad. 

Antes de introducir por el teclado los programas pueden tener necesi- 
dad de leer la sección del apéndice que describe diversos artificios de 
programación. Y finalmente, lo último pero no menos importante, ha de 
leerse el manual del usuario que se suministra con la computadora. Siem- 
pre puede encontrar algo de utilidad en su lectura. 


INGENIERIA 
ELECTRICA 


l. Parámetros de circuitos 
resonantes L-C 


DESCRIPCION DEL PROGRAMA 


Este programa evalúa un circuito resonante que esté constituido por 
una bobina de capa única, sin núcleo, sin pérdidas y un condensador. El 
circuito se describe completamente por medio de dos ecuaciones: una fór- 
mula de aproximación para la inductancia de la bobina y una relación de 
resonancia de un circuito LC. 


Fórmula 1-1 (27F) x LC = 1 _ (ND)” x 10% 


L=IEDIFAWA 


F es la frecuencia en hertzios, C es la capacidad en faradios, L es la induc- 
tancia de la bobina en henrios, A es la longitud de la bobina en pulgadas, D 
es el diámetro de la bobina en pulgadas y N es el número de espiras de la 
bobina. 

El programa ayuda a diseñar circuitos resonantes, a devanar bobinas y a 
reparar o modificar receptores de radio y transceptores. Cuando se ejecuta 
el programa, hay que introducir primero todos los parámetros conocidos y 
luego dejar que la computadora determine los que faltan por medio de las 
fórmulas anteriores. Como se muestra en los ejemplos, el programa resol- 
verá problemas bastante complicados. Los parámetros F, C, L, A,Do N 
pueden introducirse, calcularse o revisarse en orden arbitrario pulsando las 
teclas correspondientes. 
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INSTRUCCIONES 


El programa genera una petición (mensaje orientativo) de un signo de 
interrogación (?) para una orden o con el nombre de la función cuando está 
esperando un valor. Cuando haya aparecido dicho mensaje, ha de pulsarse 
una de las teclas siguientes: 


S Borrado de todas las variables, comienzo de una nueva eje- 
cución 

F Introducción de la frecuencia F en hertzios 

ZF Cálculo de la frecuencia a partir de C y de L 

O Introducción de la capacidad C en faradios 

ZC Cálculo de la capacidad a partir de F y de L 

E Introducción de la inductancia de la bobina en henrios 

LL Cálculo de la inductancia a partir de F y de C 

ZXoX Cálculo de la inductancia de la bobina a partir de A, D y N 

A Introducción de la longitud de la bobina en pulgadas 

ZA Cálculo de la longitud de la bobina a partir de L, D y N 

D Introducción del diámetro de la bobina D en pulgadas 

ZD Cálculo del diámetro de la bobina a partir de L, A y N 

N Introducción del número de espiras N para la bobina 

ZN Cálculo de números de espiras a partir de L, A y D 


Para revisar el valor de cualquier parámetros, teclear Z seguida por F, C, L, 
A,DoN. 


EJEMPLOS 


Enunciados de los problemas 


1. Diseñar un circuito resonante que utilice un condensador variable 
con una capacidad en el margen de 35 pF a 350 pF y una bobina a devanar 
en un núcleo conformador de 0,6 pulgadas de longitud y 0,3 pulgadas de 
diámetro. La frecuencia más baja debe ser de unos 6 MHz. Calcular el nú- 
mero de espiras de la bobina y la más alta frecuencia a la que resonará el 
circuito con el condensador en su valor mínimo. 


Figura 1-1 14 
espiras 


Mi ha 1.5"ao| es 
+ 
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2. Una bobina de una sola capa de 1,5 pulgadas de longitud y 0,6 pul- 
gadas de diámetro, con 14 espiras, como se indica en la figura 1-1, resuena 
con un condensador desconocido a 10 MHz. ¿Qué valor de condensador en 
paralelo se requiere para disminuir la frecuencia resonante a 9 MHz? 


Ejemplo de Ejecución 


2F ?F 
FREC.? FREC.? 
FREC.=6000000 FREC.=10000000 

20 2A 

CAP.? LONGITUD? 
CAP.=3.5E-10 LONGITUD=1.5 

27 2D 

CALCULAR DIAM. ? 

?L DIAM.=0.6 
INDUCT.=2.0103409E-6 ?N 

2A ESPIRAS? 

LONGITUD? ESPIRAS=14 
LONGITUD=0.6 2X 

2D INDUCT.=9.9661017E-7 
DIAM. ? 27 

DIAM.=0.3 CALCULAR 

E 20 

UNOS CAP.=2.5416453E-10 
: 2 
ESPIRAS=25.626381 

2N 

ESPIRAS? LET T=C 

ESPIRAS=25 

2X sE 
INDUCT.=1.9132653£-6 e 

27 : 

OD FREC.=9000000 

2F 

FREC.=6150331.3 dd 

20 j 

ALS CAP.=3.1378337E-10 
CAP.=3.5E-11 i 

27 

CALCULAR PRINT C-1 

?F 5.9618838E-11 


FREC.=19449055 
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Discusión de los resultados 


El primer problema da como resultado 25,6 espiras. En una aplicación 
práctica este número se redondea por defecto a 25 espiras. La frecuencia 
resonante es, entonces, de 6,15 MHz. La más alta frecuencia a la que el 
circuito resonará con un condensador de 35 pF es, entonces, de 19,44 
MHz. 

En el segundo problema, calculamos un condensador resonante de 254 
pF a 10 MHz. A continuación, interrumpimos la ejecución del programa 
con BREAK y almacenaremos este valor temporalmente como una varia- 
ble arbitraria T. La ejecución del programa se reanuda, entonces, con la 
tecla CONT y se determina la capacidad resonante a 9 MHz (313 pF). El 
programa vuelve a interrumpirse y la diferencia se encuentra que es de 59 
pF, como respuesta final. Este ejemplo ilustra el uso alternativo del Timex/ 
Sinclair 1000 en el modo de programación e inmediato (calculadora). 


OBSERVACIONES SOBRE LA PROGRAMACION 


El programa ilustra el empleo de la secuencia de INKEYS$ con comien- 
zo en la línea 40. La función CODE en la línea 95 proporciona un número 
(ver página 137-139 del Manual del Usuario) para cada tecla pulsada du- 
rante la secuencia de INKEYS. El flujo del programa se dirige, entonces, a 
un número de línea correspondiente a la tecla pulsada (línea 105). Si la te- 
cla no es reconocida, una rutina de intercepción de errores proporciona un 
nuevo mensaje orientativo en la línea 30. Un indicador (Z) se pone a «0» o 
a «l» dependiendo de si el programa espera una entrada (frecuencia, in- 
ductancia, etc.) o si debe encontrar un parámetro que falta a partir de los 
que se dan. 


LISTADO DEL PROGRAMA 


30 PRONTO eS 

40 IF INKEY$="" THEN GOTO 40 
45 PAUSE 10 

50 LET A$=INKEYV$ 

60 PAUSE 5 

65 PRINT AS$ 

70 IF A$=INKEY$ THEN GOTO 70 
95 LET K=CODE A$ 

100 IF K=56 THEN GOTO 20 


105 
106 


110 
380 
381 
382 
383 
384 
385 
386 
387 
400 
401 
402 
403 
404 
405 
406 
407 
410 
411 
412 
413 
414 
415 


416 
417 
430 
431 
432 
433 
434 
435 
436 
437 
490 
491 
492 
493 
494 
495 
496 
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IF K=63 OR K=40 OR K=49 OR K=43 THEN 
GOTO K*10 

IF K=41 OR K=51 OR K=38 OR K=61 THEN 
GOTO K*10 

GOTO 30 

IF Z=1 THEN GOTO 384 
PRINT "LONGITUD?" 
INPUT A 

GOTO 386 

LET Z=0 
LET-A=((D*NYTE=3)**2=18*D*L )/40/L 
PRINT "LONGITUD=";A 
GOTO 30 

IF Z=1 THEN GOTO 404 
PRINT "CAP, ?* 

INPUT C 

GOTO 406 

LET Z=0 

LET Cs TI(2APISFIESZLL 
PRINT “CAP ¿="3 

GOTO 30 

IF Z=1 THEN GOTO 414 
PRINT "DIAM.?" 

INPUT D 

GOTO 416 

LET Z=0 

LET D=(9*L+SQR (81*L*L+40*L*A*N*N/TE-6) 
) /N/N*1E6 

PRINT "DIAM.=";D 

GOTO 30 

IF Z=1 THEN GOTO 434 + 
PRANT. “FREQS25 

INPUT F 

GOTO 436 

LET Z=0 
LETF=1/2/PI/SOR AL*C) 
PRIENT:“ERECI=S "SF 

GOTO 30 

IF Z=1 THEN GOTO 494 
PRINT "INDUCT.?" 

INPUT L 

GOTO 496 

LET Z=0 

LET ¡L=IC2*P1*E)**%27c 
PRINT "INDUCT.="3;L 


5 
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497 
510 
511 
91.2 
913 
514 
915 
516 
3:17 
518 
610 
611 
612 
613 
630 
631 
632 
1000 
1010 


Originalmente publicado con el título de «La calculadora de bolsillo re- 
suelve el problema de la resonancia con el emplo de BASIC». Reimpresión 
de Electronics, 16 junio 1981, Copyright O McGraw-Hill, Inc. 1981. Todos 


GOTO 30 

IF Z=1 THEN GOTO 514 
PRINT "ESPIRAS?" 
INPUT N 

GOTO 517 

LED :Z=0 

GOSUB 1000 

LET N=SQR (X*L)/D 
PRINT "ESPIRAS=";N 
GOTO 30 

LET Z=0 

GOSUB 1000 

LEFT L=(N*D)**2Z/X 
GOTO 496 é 
LET Z=1 

PRINT "CALCULAR” 
GOTO 30 

LET X=(18*D+40*A)*1E6 
RETURN 


los derechos reservados. 


2. Redes resistivas 
de atenuación 
| y de adaptación 


DESCRIPCION DEL PROGRAMA 


Este programa calcula los valores de los componentes resistivos para re- 
des de atenuación y de adaptación que suelen utilizarse en redes de trans- 
misión. Calcula las resistencias en serie y en paralelo para seis tipos de 
redes resistivas: T y «Pi» desequilibradas, O y H equilibradas, y las redes L 
desequilibrada o C equilibrada adaptadoras de pérdidas mínimas. Además 
de calcular los valores de las resistencias, el programa determina también 
las pérdidas mínimas de una red adaptadora de impedancias. Si no se re- 
quiere un aislamiento de corriente continua, una red adaptadora de impe- 
dancias suele utilizarse como un sustituto sencillo y barato de un transfor- 
mador adaptador de impedancias. La designación de los componentes se 
muestra en las figuras, 2-1, 2-2, y 2-3. Si el programa determina valores ne- 
gativos para las resistencias para una impedancia de entrada, una impedan- 


Figura 2-1 Red T desequilibrada/H equilibrada RI/2 R2/2 


Ri r2 
R3 
Z1| — —— 22 
R3 
RI/2 rR2/2 
N K 1 () PERDIDAS (dBy 10 R2 = Z2 N_+ ] = R3 
NN *=" 1 
2 VNZI x Z2 
E — BERET O 
N=1 RI = ZI (7) — R3 
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REDES DE ATENUACION 

T DESEQ., PI EQUIL. - INTROD.1 
PI DESEQ., O EQUIL. - INTROD.2 
PERDIDAS MIN. - INTRODUCIR 3?2 
INTRODUCIR Z1?300 

INTRODUCIR Z2?300 

INTRODUCIR PERDIDAS EN DB?10 
R1=577.48518 OHM 

R2=577.48518 0HM 

R3=426.90748 0HM 


REDES DE ATENUACION 

T DESEQ., PI EQUIL. - INTROD.] 
PI DESEQ., O EQUIL. -.INTROD.2 
PERDIDAS MIN. - INTRODUCIR 32?3 
TVS" 7e 

INTRODUCIR Z1?300 

INTRODUCIR Z2?75 
PERDIDAS=0.57194755 DB 
R1=259.80762 OHM 

R3=86.60254 0OHM 


REDES DE ATENUACION 

T DESEQ., PI EQUIL. - INTROD.1 
PI DESEQ., O EQUIL. - INTROD.2 
PERDIDAS MIN. - INTRODUCIR 3? 


Discusión de los resultados 


Por regla general, los valores de resistencias normalizados pueden en- 
contrarse dentro de una pequeña tolerancia porcentual de los valores 
calculados. Obsérvese que, en el tercer ejemplo, las pérdidas mínimas son 
de sólo 0,57 dB. Ello podría compararse con el valor entre 1 y 2 dB para un 
transformador de adaptación. 


OBSERVACIONES SOBRE LA PROGRAMACION 


Las comprobaciones de errores para los valores de entrada se realizan 
en la línea 75. A continuación, las fórmulas adecuadas son objeto de selec- 
ción dependiendo del tipo de red elegido. Para conservar el código se utili- 
za la misma secuencia de entrada (líneas 90-125) para Z1 y Z2 para los tres 
tipos de redes. Análogamente, la misma secuencia de salida (impresión) (lí- 
neas 350-370) para las tres ramas. 


REDES RESISTIVAS DE ATENUACION Y DE ADAPTACION 


LISTADO DEL PROGRAMA 

10. CES 

30 PRINT "REDES DE ATENUACION" 
40-“PRINT.T-DESEOS<5/ PT EQUILS +=. INTROD: 1" 
50 PRINT.*PT DESEO... 0 EQUIL. y INTROD. 2* 
60 PRINT "PERDIDAS MIN. - INTRODUCIR 32"; 
70 INPUT A 

72 PRINT A 

75 IF A<>INT A OR A>3 OR A<] THEN GOTO 280 
80 IF A=3 THEN PRINT "Z1l > Z2" 

90 PRINT "INTRODUCIR Z1?"; 

100 INPUT B 

105 PRINT B 

110 PRINT "INTRODUCIR Z2?"; 

120 INPUT C 

125 PRINT C 

130 IF A=3 THEN GOTO 270 

140 PRINT "INTRODUCIR PERDIDAS EN DB?"; 
150 INPUT D 

155 PRINT D 


LET N=10**(D/10) 

IF A>1 THEN GOTO 220 
LET G=2*SQR (N*B*C)/( 
LET E=B*(N+1)/(N-1)-G 
LET F=C*(N+1)/(N-1)-G 
GOTO 350 

IF A=3 THEN GOTO 270 
LET G=(N-1)/2*SQR (B*C/N) 

LET E=1/((N+1)/(N-1)/B-1/G) 

LET F=1/((N+1)/(N-1)/C-1/G) 

GOTO 350 

IF 1-C/B>0 THEN GOTO 300 

PRINT "ERROR, REINTRODUCIR" 

GOTO 40 

LET E=B*SQR (1-C/B) 

LET G=C/SQR (1-C/B) | 

LET L=LN (SQR (B/C)+SQR (B/C-1))/LN 10 
PRINT "PERDIDAS=";L;" DB" 

PRINT "Rl=";E;" OHM" 

IF A<>3 THEN PRINT "R2=";F;" OHM" 
PRINT "R3="3G;" OHM" 

PRINT ** 

GOTO 20 


N-1) 
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3. Impedancia característica 
de las lineas de transmisión 


DESCRIPCION DEL PROGRAMA 


Este programa calcula la impedancia característica correspondiente a 
los tipos más comunes de líneas de transmisión: el conductor coaxial sim- 
ple, la línea blindada equilibrada, la línea bifilar abierta en aire, la línea de 
bandas paralelas utilizada a las frecuencias de microondas y la línea de 
transmisión constituida por dos' hilos blindados en paralelo con retorno 
provisto de vaina. Los cinco tipos de líneas de transmisión se muestran en 
las cinco figuras siguientes con las fórmulas correspondientes. Las unidades 
de longitud son arbitrarias (pulgadas, pies, centímetros, etc.), pues sólo re- 
laciones de medidas lineales se utilizan en las fórmulas de cálculo. 


l. Línea coaxial simple 


Figura 3-1 
138 D 
ZO = log |5 
| Ve > (2 
R D e: constante dieléctrica relativa 


de 


2. Línea blindada equilibrada 


Figura 3-2 
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3. Línea bifilar abierta en aire 


Figura 3-3 


> ollo VR 
a 


4. Línea de bandas paralelas 


Figura 3-4 
E ii Y ri W 

€ 
ls al W | H/W < 0.1 


5. Dos hilos en paralelo con retorno provisto de vaina 


Figura 3-5 La misma figura que figura 3-2 
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LISTADO DEL PROGRAMA 
15" CES 
20 PRINT "IMPEDANCIA CARACT." 
25 PRINT 


27 PRINT "INTRODUCIR TIPO (1-5)?"; 

30 INPUT A 

35 PRINT A 

37 1F A<>INT A OR A<] OR A>5 THEN GOTO 60 
40 PRINT “INTROD.*CONST.. DIEE:- REL3?": 
45 INPUT E 

50 PRINT E 

55 IF E>0 THEN GOTO 100%*A 

60 LET L=25 

65 GOTO 720 

100 PRINT "LINEA COAXIAL SIMPLE" 

110 GOSUB 650 

120 IF R>0 AND D>R THEN GOTO 150 

130 LET L=110 

140 GOTO 720 

150 LET Z=138/SQR E*LN (D/R)/LN 10 

160 GOTO 580 

200 PRINT "LINEA BLINDADA EQUILIBRADA" 
210 GOSUB 600 

220 IF R>0 AND H>R AND D>HA+R THEN GOTO 250 
¿30. LET L=ZH0 

240 GOTO 720 

250 LET V=H/R 

260 LET S=H/D 

2JO.-LET Z>276S0R EEN LANAS 

(1+5*S))/LN 10 

280 GOTO 580 

300 PRINT "LINEA BIFILAR” 

310 GOSUB 600 

320 IF R>0 AND H>R THEN GOTO 350 

330 LET .L=310 

340 GOTO 720 

350 LET S=H/R 

360 LET Z=276*LN (S+SQR (S*S-1))/LN 10/SQR E 
3/0 GOTO 580 

400 PRINT "LINEA BANDAS PARALELAS" 
410 PRINT "ANCHURA CONDUCTOR W?"; 
420 INPUT QY 

430 PRINT Y 
440 GOSUB 600 
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450 IF>H O AND W>0 THEN GOTO 480 

460 LET L=410 

470 GOTO 720 

480 LET Z=377*H/W/SQR E 

490 GOTO 580 

500 PRINT "2 LINEAS PAR., RET. CON VAINA" 
510 GOSUB 600 

520 IF R>0 AND H>R AND D>H+R THEN GOTO 550 
930 LET L=510 

540 GOTO 720 

550 LET V=H/R 

560 LET S=H/D 

SO LET 2=69/SQR. ELN (V23S* (TES FFAJOSUEN TO 
580 PRINT: “ZO0="52Z3%710HM" 

590 GOTO 25 

600 PRINT "SEPARACION H?"; 

610 INPUT H 

620 PRINT H 

630 IF A=3 THEN GOTO 680 

640 IF A=4 THEN RETURN 

650 PRINT "DIAM. EXT. D?"; 

660 INPUT D 

670 PRINT D 

680 PRINT "DIAM. COND. R?"; 

690 INPUT R 

700 PRINT R 

710 RETURN 

720 PRINT "ERROR, REINTRODUCIR" 

730 GOTO L 


TRANSMISION 
DE DATOS 


4. Tasas de errores 
en bits para diversos 
planes de modulación 


DESCRIPCION DEL PROGRAMA 


En un enlace de transmisión de datos bien diseñado, el ruido aleatorio 
con distribución gausiana es la fuente principal de errores de bits de datos. 
La probabilidad de recibir un bit de datos incorrecto depende de la relación 
señal/ruido (S/N) y del tipo de modulación de datos utilizada. Por regla ge- 
neral, la complejidad del circuito está ligada con la probabilidad de error. 
Un circuito más complejo, con más estrictas exigencias de atenuación y de 
distorsión de retardo dará lugar a una menor probabilidad de error para la 
misma relación S/N. El siguiente programa calcula la probabilidad de error 
de bit (Pe) como una función de S/N para nueve planes de modulación 
bien cómodos. Las fórmulas utilizadas son las siguientes: 


1. Manipulación todo-nada (On-Off Keying = OOK), coherente: 


Fórmula 4-1 Pe = 1 Erfe (Y y) 


2. Manipulación por desplazamiento de frecuencia (FSK), coherente 
Manipulación por desplazamiento de amplitud (ASK), coherente 
Modulación por codificación de impulsos (PCM), unipolar 


Aa Ss 
Erfc 2 ÑN 


Fórmula 4-2 Per= 


n|i— 
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3. PCM, polar 4. FSK, no coherente 
Manipulación por desplazamiento de 
fase (PSK) 
Fórmula 4-3 Fórmula 4-4 
e. El -1.(5%) 
Pe = > Erfc N Pe =38€ 


5. Manipulación por desplazamiento de 6. ASK, no coherente 
fase diferencial (DPSK) 


Fórmula 4-5 AR 


Pe =>3el-%) pe = Emo [2Y y) 


La función Erfc (función de error complementaria) relacionada con la distribución 
de ruido gaussiano se define como sigue. La función se calcula de forma recurrente, 
como se indica en el programa 12. 


Fórmula 4-7 Erf (x) = Ly er” de Erfc (x) = 1 — Erf (x) 
0 


INSTRUCCIONES 


Cuando se pulsa la tecla RUN, el programa mide la relación señal/rui- 
do en dB y la clase de modulación. Seleccionar una relación S/N y un 
número entre 1 y 6 para indicar la clase de modulación tal como se descri- 
be en la descripción del programa. Entonces, el programa calcula la 
probabilidad de error de bits. 


EJEMPLOS 


Enunciado de problema 


1. Calcular la probabilidad de error de bit Pe para una S/N de 12 dB 
para nueve planes de modulación. 
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Ejemplo de ejecución 


PROBABILIDAD DE ERROR DE BIT 

INTROD. RELACION S/N EN DB?12 
SELECCIONAR CODIGO MODUL.1-6?1 
PROB(12 DB)=.0024385373 


INTROD. RELACION S/N EN DB?12 
SELECCIONAR CODIGO MODUL.1-6?2 
PROB(12 DB)=.000034302624 


INTROD. RELACION S/N EN DB?12 
SELECCIONAR CODIGO MODUL.1-6?3 
PROB(12 DB)=9.0060101E-9 


INTROD. RELACION S/N EN DB?12 
SELECCIONAR CODIGO MODUL.1-67?4 
PROB(12 DB)=.0001808915 


INTROD. RELACION S/N EN DB?12 
SELECCIONAR CODIGO MODUL.1-6?5 
PROB(12 DB)=6.544347E-8 


INTROD. RELACION S/N EN DB?12 


SELECCIONAR CODIGO MODUL.1-6?6 
PROB(12 DB)=.0001990186 


Discusión de resultados 


La probabilidad de error para los nueve planes de modulación varía 
mucho. La más baja probabilidad de error se da por el plan de modulación 
tipo 3 (PCM, Polar/PSK) y la más alta por el plan de modulación tipo 1 


(OOK, coherente). 


OBSERVACIONES SOBRE LA PROGRAMACION 


Las fórmulas aparentemente diferentes se resuelven por la misma sec- 
ción de código recalculando el argumento de la función de error comple- 
mentaria. Una nueva discusión de evaluación recurrente de la función de 


error se da en el programa 12. 


LISTADO DEL PROGRAMA 
15 CLS 


20 PRINT "PROBABILIDAD DE ERROR BIT" 
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PRINT "INTROD. RELACION S/N EN DB?"; 


INPUT B 
PRINT B 
LET C=B 


PRINT "SELECCIONAR CODIGO MODUL.1-6?*; 


INPUT D 
PRINT D 


IF D>0 AND D<7 AND D=INT D THEN GOTO 100%*D 


PRINT "ERROR ENTRADA" 
GOTO 45 

GOSUB 800 

GOSUB 800 
LET-A=LO**1C720) 

LET Y=1/SQR PI*EXP (-X*X) 
IF X<2 THEN GOTO 350 
LET A=14/X 

FOR I=27 TO 1 STEP -1 
LET A=1/2/(X+A) 

NEXT 1 

LET Z=Y/(X+A) 

LET N=1-Z 

GOTO 700 

LET “11 

LET T=2*X*Y 

LET N=T 

LET T=2*X*X*T/(2*1+1) 
LET M=N+T 

IF M=N THEN GOTO 393 
LET N=M 

LET 1=1+1 

GOTO 365 

LET Z=1-N 

GOTO 700 

GOSUB 800 

LET-ZSEXP =(T10**(€/103) 
GOTO 700 

LEFT. F=10**1C7/10:) 


LETOLSERP LFL2IALIELISOR EEAPIAF 


LET Z=ZP2 

PRINT "PROB(";B;" DB)=";Z 
PRINT 

GOTO 25 

LET C=C-10*LN 2/LN 10 
RETURN 


TEORIA 
DE LOS NUMEROS 


5. Conversión de cualquier 
base a cualquier otra base 


DESCRIPCION DEL PROGRAMA 


La conversión de números entre diferentes bases se suele realizar por 
programadores de computadoras, matemáticos, estudiantes de segunda en- 
señanza y universitarios. En particular, el desarrollo de un programa de 
computadora exige la conversión entre las bases de numeración binaria (2), 
octal (8), decimal (10) y hexadecimal (16). Los problemas de la teoría de 
números pueden requerir la conversión entre éstos y otras bases. 

La conversión entre dos bases arbitrarias a y b suele ser un procedi- 
miento de dos pasos, que convierte primero desde la base a a la base 10 y 
luego, desde la base 10 a la base b. El método singular aquí presentado para 
enteros positivos utiliza el mismo algoritmo para ambas conversiones, con 
el ahorro consiguiente en pasos de programación. 

Para representar un número en una base más grande que 10, se necesita 
más de una posición para representar un dígito en esa base. Por ejemplo, 
17 FE (base 16) se introduciría o visualizaría como [01], [07], [15], [14] o 
simplemente 1071514 (1 =01, 2=02,..., A=10, B=11, C=12, D=13, 
E=14, F=15). Por lo general, una sola posición se reservaría por dígito 
para las bases 2-10, dos posiciones para bases 11-100, etc. El programa in- 
terpreta, de forma automática, la entrada o la visualización dependiendo 
de la base. El programa no tiene ninguna limitación inherente sobre la 
magnitud de la base o el número a convertir. El único factor limitador es la 
precisión de la computadora. Por ejemplo, la precisión del Timex/Sinclair 
BASIC limitará la entrada y la salida a unas 10 posiciones decimales. Un 
número más grande daría lugar a un error. 


INSTRUCCIONES 


Después de pulsar la tecla RUN, hay que atender las peticiones del pro- 
grama introduciendo el número a convertir (argumento) y las bases antigua 
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y nueva. El programa visualizará el argumento convertido a la nueva base; 
si ni la base antigua, ni la nueva, es igual a 10, entonces, se visualizará tam- 
bién el argumento en la base 10. 


EJEMPLOS 


1. Convertir 17C (base 16) a base 2. 

2. Convertir el equivalente decimal de 17C a base 2 (ambos resultados 
deben dar la misma respuesta). 

3. Convertir el argumento decimal 12345 a base 122. 

4. Convertirlo de nuevo desde base 122 a base 10. 


Ejemplo de ejecución 


CONVERSION BASE GENERAL 
INTROD. ARGUMENTO?10712 
INTROD. BASE ANTIGUA?16 
INTROD. BASE NUEVA?2 
10712 A BASE 16 
=1011111 A BASE 2 

=380 A BASE 10 


INTROD. ARGUMENTO?380 
INTROD. BASE ANTIGUA?10 
INTROD. BASE NUEVA?2 
380 A BASE 10 
=101111100 A BASE 2 


INTROD. ARGUMENTO? 12345 
INTROD. BASE ANTIGUA?10 
INTROD. BASE NUEVA?122 
12345 A BASE 10 
=101023 A BASE 122 


INTROD. ARGUMENTO?101023 
INTROD. BASE ANTIGUA?122 
INTROD. BASE NUEVA?16 
101023 A BASE 122 
=3000309 A BASE 16 
=12345 A BASE 10 
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Discusión de los resultados 


Se requiere alguna interpretación de los resultados para bases mayores 
que 10. El resultado del tercer ejemplo debe interpretarse como [101], 


[023] (base 122). 


OBSERVACIONES SOBRE LA PROGRAMACION 


Una amplia verificación de errores de los datos de entrada se realiza en 
las líneas, 70, 110 y 150. El error de redondeo se hace mínimo utilizando la 
función INT y sumando 0,5 a las funciones adecuadas en las líneas 440 y 
480. La línea 480 determina el número de dígitos decimales requeridos 
para representar cada «dígito» para una base arbitraria. 


LISTADO DEL PROGRAMA 

20 CLS 

30 PRINT "CONVERSION BASE GENERAL” 
40 PRINT "INTROD. ARGUMENTO?"; 

50 INPUT C 

60 PRINT C 

70 IF C<>INT (ABS C) THEN GOTO 40 
80 PRINT "INTROD. BASE ANTIGUA?"; 
90 INPUT D 

100 PRINT D 


IF D<>INT (ABS D) THEN GOTO 80 
PRINT "INTROD. BASE NUEVA?"; 


INPUT E 
PRINT E 


IF E<>INT (ABS E) THEN GOTO 120 


IF D<>10 THEN GOTO 190 
LET N=C 

GOTO 250 

LET U=D 

GOSUB 480 

LET S=C 

LET Q=D 

LET R=V 

GOSUB 390 

IF E<>10 THEN GOTO 280 
LET S=N 

GOTO 340 

LET U=E 

GOSUB 480 


300 
310 
320 
330 
340 
390 
360 


370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
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LET S=N 

LET Q=V 

LET R=E 

GOSUB 390 

PRINTS" A: BASED 

PRINT "=";3N;5" A BASE ";E 

IF D<>]0 AND E<>10 THEN PRINT *=";S;" 
TO BASE 10" 

PRINT + 

GOTO 40 

LET 
LE-T 
CET 
LET 
LET NT (P/R) 

LET N=N+INT ((T-P*R)*Q**M+.5) 

IF P=0 THEN RETURN 

LET M=M+] 

GOTO 420 

LET V=INT (10**(1+INT (LN (U-1)/LN 10))+.5) 
RETURN 


UU AUZZ 
A A 


HH TUDONOO 


6. Conversión hexadecimal 
a decimal y viceversa 


DESCRIPCION DEL PROGRAMA 


Este programa difiere del anterior programa de conversión de cualquier 
base en que sólo convierte entre bases 10 y 16. Las dos conversiones (10 a 
16 y 16 a 10) son muy importantes para desarrollo de programa de compu- 
tadora, pues la mayor parte de los dispositivos de programación (ensambla- 
dores, mapas de memoria, desensambladores, etc.) se refieren a las posicio- 
nes de memoria y al contenido de memoria en base 10 ó 16 (hex). Este 
programa, a diferencia con el anterior, utiliza la presentación alfanumérica 
estándar para dígitos hexadecimales: O a 9 inclusive y A, B, C, D, E y F, 
para los dígitos 10 a 15 inclusive. Sigue siendo aplicable el comentario 
acerca de la precisión y del redondeo del programa anterior. 


INSTRUCCIONES 


Pulsar la tecla RUN y seguir las indicaciones (mensajes orientativos) 
para realizar la conversión hexadecimal-decimal (1) o decimal-hexade- 
cimal (2). Después de la siguiente indicación, introducir el argumento en 
hexadecimal (empleando A a F, si fuera necesario) o como un número de- 
cimal. 


EJEMPLOS 


1. Convertir F697 (base 16) a base 10. 
2. Convertir el resultado de nuevo a base 16. 
3. Convertir FFFF (base 16) a base 10. 


Ejemplo de ejecución 


HEX A DECIMAL, INTROD.1 
DECIMAL A HEX, INTROD.2?1 
INTROD. ARGUMENTO EN HEX?F697 
=63127 EN BASE 10 


HEX A DECIMAL, INTROD.1 
DECIMAL A HEX, INTROD.2?2 
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CONVERSION HEXADECIMAL A DECIMAL Y VICEVERSA 


INTROD. ARGUMENTO A BASE 10263127 
=F697 EN BASE 16 


HEX A DECIMAL, INTROD.1 
DECIMAL A HEX, INTROD.2?1 
INTROD. ARGUMENTO EN HEX?FFFF 
=65535 EN BASE 10 


Discusión de los resultados 
Las indicaciones y los resultados son auto-explicatorios. 


OBSERVACIONES SOBRE LA PROGRAMACION 
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La comprobación de errores de los datos de entrada se realiza en las lí- 
neas 70 y 270. Los dígitos en base 16 se introducen como una cadena de 
caracteres (A$) y cada carácter en la cadena se interpreta en la línea 260. Si 
se encuentra un carácter que no corresponde a un número entre O y 15, se 
rechaza como un error. La traducción desde números decimales a una ca- 
dena de caracteres en hexadecimal se produce en la línea 510. Este progra- 
ma pone de manifiesto la potencia de las funciones de manipulación de 


cadenas de caracteres en BASIC, tales como CHR$, CODE y LEN. 


LISTADO DEL PROGRAMA 


10 DIM B$(10) 

20 CLS 

30 PRINT | 

35 PRINT "HEX A DECIMAL, INTROD.1" 
40 PRINT "DECIMAL A HEX, INTROD.2?"; 
50 INPUT B 

60 PRINT B 

70 IF B<>I AND B<>2 THEN GOTO 30 

80 GOTO 200%*B 

200 PRINT "INTROD. ARGUMENTO EN HEX2?"; 
210 INPUT AS 
220 PRINT AS 
230 LET L=LEN AS 

240 LET N=0 
250 FOR 1=0 TO L-1 
260 LET C=CODE A$(L-1)-28 

270 IF C>=0 AND C<=16 THEN GOTO 310 
280 PRINT "ENTRADA ERRONEA" 
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290 
300 
310 
320 
330 
340 
400 
410 
420 
430 
440 
450 
460 
470 
480 
485 
490 
500 
510 
580 
590 
600 
610 
620 
630 
640 
650 
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PAUSE 200 

RUN 

LET N=N+C>16**1 

NEXT 1 | 

PRINT. *S"SNs” EN. BASE “10 

GOTO 640 

PRINT "INTROD. ARGUMENTO A BASE 102"; 


INPUT D 
PRINT D 


IF D=INT (ABS D) THEN GOTO 470 
PRINT "ENTRADA ERRONEA” 

PAUSE 200 

RUN 

LET K=0 

LET K=K+1 

LET F=D 

LET D=INT (D/16) 

LET E=F=.16*B 

LET B$(K)=CHR$ (E+28) 


IF F>0 THEN GOTO 480 
PRINT "="; 

FOR I=1 TO K-1 

PRINT BS$(K-I1); 

NEXT 1 

PRINT "EN BASE 16" 
PRINT " " 

GOTO 30 


7. Descomposición 
en factores primos 


DESCRIPCION DEL PROGRAMA 


Uno de los primeros problemas con los que se enfrenta un estudiante en 
la teoría de los números es la descomposición en primos de un entero. Esta 
Operación no es simplemente un ejercicio teórico (es aplicable en muchas 
áreas incluyendo la criptografía). Algunos de los modernos códigos cripto- 
gráficos están basados en el supuesto de que lleva mucho tiempo, incluso 
para computadoras grandes, determinar los factores primos de números 
grandes. En la técnica denominada «clave pública», puede publicarse el 
código de cifrado (producto de dos primos), mientras que los factores pri- 
mos que constituyen el código de descifrado son conocidos solamente para 
los receptores autorizados. Si un código está basado en una clave que es un 
producto de dos primos grandes, lo que llevaría años a una computadora 
determinar, el código será inviable a todos los efectos prácticos. 

Los dos métodos principales para encontrar los factores primos de un 
entero son la criba de Eratóstenes, así llamada en honor a un famoso mate- 
mático griego y una sencilla comprobación de divisores impares. El primer 
método es más rápido pero también más adecuado para grandes computa- 
doras, pues almacena en memoria todos los primos sucesivos más peque- 
ños que la raíz cuadrada del número a descomponer. El segundo método es 
mucho más sencillo, pues solamente comprueba la descomposición en 2, 
3, 5 y los sucesivos números enteros impares. Aunque el segundo método 
desperdicia algo de tiempo de ejecución, puesto que comprueba todos los 
números impares y no sólo los números enteros primos, tiene menos nece- 
sidades de memoria en la computadora, porque no han de almacenarse 
resultados intermedios. Por consiguiente, el segundo método es el único 
adecuado para el cálculo de números primos en la computadora Timex/ 
Sinclair 1000 y es la base para el programa siguiente. 


INSTRUCCIONES 


Pulsar la tecla RUN para iniciar el programa y luego, seguir las indica- 
ciones introduciendo el número a descomponer en factores primos. A 
continuación, el programa visualiza cuantos números primos se encontra- 
ron y dichos factores primos, uno por uno. 
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EJEMPLOS 


Enunciado del problema 


Descomponer 12345, 66759 y 71 en factores primos. 
Ejemplo de ejecución 


BUSCADOR FACTORES PRIMOS 
INTRODUCIR NUMERO?12345 
NO. DE FACTORES=3 

FACTOR NO. 1=3 

FACTOR NO. 2=5 

FACTOR NO. 3=823 


BUSCADOR FACTORES PRIMOS 
INTRODUCIR NUMERO?66759 
NO. DE FACTORES=5 


FACTOR NO. 1=3 
FACTOR NO. 2=7 
FACTOR NO. 3=11 
FACTOR NO. 4=17 
FACTOR NO. 5=17 


BUSCADOR FACTORES PRIMOS 
INTRODUCIR NUMERO?71 

NO. DE FACTORES=1 
FACTOR NO. 1=71 


Discusión de resultados 


Obsérvese, en el segundo ejemplo, que 17 aparece dos veces, lo que in- 
dica que es un factor múltiple. En el tercer ejemplo, probamos un número 
primo (71) que tiene un factor no trivial único, que es el mismo. 


OBSERVACIONES SOBRE PROGRAMACION 


El programa realiza la comprobación de errores en la línea 90. Si el nú- 
mero a descomponer es más grande que 1.000.000.000, los errores del 
redondeo no permitirán que el programa actúe correctamente. Para llama- 
da del usuario los cincuenta primeros factores primos están almacenados 
en una matriz («array») dimensionada B(50) y podrían llamarse por medio 
de un programa escrito por el usuario. Si hay más de 50 factores, el progra- 


ma visualizará solamente los cincuenta primeros. 
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LISTADO DEL PROGRAMA 


20; €L3 

25 LET K=0 

27 LET D=0 

30 DIM B (50) 

40 PRINT "BUSCADOR FACTORES PRIMOS" 
50 PRINT "INTRODUCIR NUMERO?”"; 

60 INPUT A 

70 PRINT A 

80 LET F=A 

90 IF A=INT A AND A>1 AND A<1E9 THEN.GOTO 120 
100 PRINT "ERROR, REINTRODUCIR" 

110 GOTO 50 

120 LET E=A/2 

130 LET C=2 

140 IF E<>INT E OR 2*E<A THEN GOTO 170 
150 GOSUB 320 

160 GOTO 120 

170 LET D=D+1 

180 LET C=1+2*D 

190 IF C>SQR A THEN GOTO 250 

200 LET E=A/C 
210 IF E<>INT E OR E*C<>A THEN GOTO 170 
220 LET D=D-1 
230 GOSUB 320 

240 GOTO 170 
250 IF A<=1 THEN GOTO 280 
260 LET K=K+1 
270 LET B(K)=A S 
280 PRINT " NO. DE FACTORES=";K 
290 FOR I=1 TO K 

300 PRINT: “EACTOR NO. "31:35" B1 1) 
310 NEXT 1 
312 PRINT "2 

IS “GOTUO.: 23 
320 IF K>=50 THEN RETURN 
330 LET K=K+1 
340 LET B(K)=C 
350 LET A=E 
360 RETURN 


PROGRAMACION 
DE COMPUTADORA 


8. Cálculos del 
desplazamiento (offset) 
relativo para programación 


en código máquina 


DESCRIPCION DEL PROGRAMA 


La finalidad de este programa es calcular las direcciones relativas de 
instrucciones en código máquina. Cuando se escriben programas cortos en 
código máquina para cualquiera de los microprocesadores de 8 bits más co- 
nocidos, tales como Z80, 6502 o 6800, uno de ellos ha de cargar el progra- 
ma ensamblador, que puede no estar disponible, o ensamblar «a mano» el 
programa en código máquina. El problema principal es, entonces, encon- 
trar los argumentos de las instrucciones de bifurcación relativa. El progra- 
ma aquí mostrado determinará, en notación hexadecimal, el argumento de 
una instrucción de bifurcación relativa ZZ dada la dirección absoluta de 
esa instrucción, la dirección «FROM» (X) y la dirección hexadecimal de la 
instrucción a la que el programa debe bifurcar si es verdadera la condición 
de bifurcación, la dirección «TO» (Y). Por ejemplo, dado el siguiente códi- 
go en lenguaje máquina del Z80: 


Dirección Instrucción Rótulo Nemónico 
20FGC DD7E 00 LOOP LD A,(1X) 
2109 B8 CPB 
210A 20 zz JRNZLOOP 
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se ha de encontrar el valor de zz, también conocido como desplazamiento 
(offset), dada la dirección de la instrucción de bifurcación «FROM» 210A 
(hex) y la dirección de la instrucción «TO»20FC (hex). Para calcular zz, 
debe resolverse una de las ecuaciones siguientes: 


zz= Y - X -2 para Y > X (bifurcación hacia adelante) o 
zz = FF (hex) + Y - X-— 1 para Y < =X (bifurcación hacia atrás) 


El valor offset resultante en el anterior ejemplo sería FO (hex) y la últi- 
ma instrucción leería 20F0. Habida cuenta de que el offset zz sólo puede 
tomar valores entre O y FF (hex) para microprocesadores de 8 bits, la bifur- 
cación relativa está limitada a media página (1 página = 256 posiciones) en 
el sentido hacia adelante o hacia atrás. El programa comprueba si se puede 
realizar la instrucción de bifurcación, pues analiza las direcciones de las 
instrucciones FROM/TO. El programa indicará si la bifurcación relativa 
puede, o no puede, realizarse en un solo paso y si dará lugar a un salto ha- 
cia adelante o hacia atrás. 


INSTRUCCIONES 
Cuando se pulsa la tecla RUN, el programa realiza la llamada de la di- 


rección de las instrucciones «FROM» y «TO». Ambas deben introducirse 
en notación hexadecimal. 


EJEMPLOS 


Enunciado del problema 


Instrucción ' Instrucción de 
bifurcación dirección destino dirección 
Ejemplo n.? absoluta (hex) absoluta (hex) 
1 123A 127B 
2 25F3 2613 
3 010A 00FC 
4 AA25 BA24 


Ejemplo de ejecución 
DIRECCIONAMIENTO REL. (8 BITS) 


INTROD. DIREC. "FROM" EN HEX?123A 
INTROD. DIREC. "TO" EN HEX?1278B 
OFFSET ADELANTE=3F 
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INTROD. DIREC. "FROM" EN HEX?25F3 
INTROD. DIREC. "TO" EN HEX?2613 
OFFSET ADELANTE=1E 


INTROD. DIREC. "FROM" EN HEX?10A 
INTROD. DIREC. "TO" EN HEX?FC 
OFFSET ATRAS O EN LUGAR=FO 


INTROD. DIREC. "FROM" EN HEX?AA25 
INTROD. DIREC. "TO" EN HEX?BA24 
4095(DEC) OFFSET DEM. GRANDE,REINT. 


INTROD. DIREC. "FROM" EN HEX?AA25 
INTROD. DIREC. "TO" EN HEX?AA23 
OFFSET ATRAS O EN LUGAR=FC 


Discusión de resultados 


Los tres primeros ejemplos dan lugar a offsets hexadecimales de 3F, 
1E y FO. El cuarto ejemplo indica que el offset es demasiado grande. Nos 
percatamos, entonces, de que la dirección «TO» debe ser AA23 en lugar de 
BA24 y finalizar con un offset de FC. 


OBSERVACIONES SOBRE PROGRAMACION 


Obsérvese el empleo de dobles comillas (Q desplazada) en las líneas 45 
y 120. Para realizar cálculos aritméticos en número hexadecimales, los nú- 
meros se convierten primero a la base 10 y los resultados se vuelven a 
convertir a base 16. La conversión de decimal a hexadecimal, y viceversa, 
se realiza de la misma forma que en el programa n.* 6 por medio de funcio- 
nes de manipulación de cadenas. 


LISTADO DEL PROGRAMA 


¿0 ELS 

30 PRINT "DIRECCIONAMIENTO REL. (8 BITS)" 
40 PRINT 

45 PRINT "INTROD. DIREC. ""FROM"" EN HEX?"; 
50 INPUT C$ 

60 PRINT C$ 

70 GOSUB 330 

80 IF F=0 THEN GOTO 110 


90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
245 
250 
260 
265 
270 
280 
290 
300 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
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PRINT "ERROR, REINTRODUCIR" 

GOTO 40 

LET J=X 

PRINT "INTROD. DIREC. ""TO"" EN HEX?"; 
INPUT C$ 

PRINT C$ 

GOSUB 330 

IF F=0 THEN GOTO 190 

PRINT "ERROR, REINTRODUCIR" 

GOTO 120 

LET D=X-J 

IF D<=129 AND D>=-126 THEN GOTO 230 


PRINT D;"(DEC) OFFSET DEM. GRANDE,REINT. 


GOTO 40 
IF D<2 THEN GOTO 260 

LET Z=D-2 

LET H$="FORWARD" 

GOTO 270 

LET Z=254+D 

LET H$="ATRAS O EN LUGAR" 
LET G=INT (Z/16) 

LET E$=CHR$ (G+28) 

LET W$=CHR$ (Z-16*G+28) 
PRINT "OFFSET "; H$;"=";ES; WS 
GOTO 40 

LEF F=0 

LET X=0 

LET L=LEN C$ 

FOR I=0 TO L-1 

LET R=CODE C$(L-1)-28 

IF R<0 OR R>15 THEN LET F=] 
LET X=X+R*16**] 

NEXT 1 

RETURN 
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TRAZADOS 
GENERADOS 
POR COMPUTADORA 


9. Alisado de curvas e 
interpolación/extrapolación 
con polinomios 
de Lagrange 


DESCRIPCION DEL PROGRAMA 


Este programa es el equivalente de una «curva francesa». Su aplicación 
principal es la interpolación y extrapolación de resultados medidos, dadas 
X(1), YG) (1 =1...., n), el programa determina un polinomio único de grado 
(n — 1) que se ajusta a estos n puntos y evalúa el polinomio para valores ar- 
bitrarios de la variable independiente X. El método difiere de la aproxima- 
ción de los mínimos cuadrados en que el polinomio resultante se ajusta 
exactamente a los puntos X(1), Y (1). 

La siguiente fórmula, atribuida al matemático francés Lagrange, descri- 
be el polinomio 


Fórmula 9-1 Te) 
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Para aplicaciones prácticas es recomendable tener cuidado en seleccio- 
nar solamente puntos contiguos para calcular el polinomio de modo que 
no se obtengan curvas con bruscas oscilaciones como extrapolaciones e in- 
terpolaciones de funciones relativamente «suaves». Por regla general, un 
polinomio de tercer grado, basado en cuatro puntos cercanos (n = 4), será 
satisfactorio. 

Una característica importante del programa es que una «ventana» de 
interpolación móvil puede establecerse de modo que sólo k valores medios 
(k < =n) más próximos al punto de interpolación se incluirán en el cálculo. 


INSTRUCCIONES 


Al pulsar la tecla RUN se inicializa el programa y emite llamadas 
orientativas correspondientes al número de puntos (n < =20) y los sucesi- 
vos valores X(1) e Y(1) de puntos en orden creciente de X, el número de 
puntos en la «ventana» y el valor de X en el que la función debe interpo- 
larse o extrapolarse. La siguiente opción es recalcular la función para otra 
X, cambiar la anchura de la «ventana» o iniciar una nueva ejecución. 


EJEMPLO 


Enunciado del problema 


Como se muestra en la figura 9-1, se han medido seis puntos. Suponer 
una anchura de «ventana» de cuatro y evaluar la función en X =2.1 (inter- 
polación) y en X =7.0 (extrapolación). A continuación, cambiar la anchu- 
ra de «ventana» a cinco puntos y volver a calcular los mismos dos puntos. 


Ejemplo de ejecución 


INTERPOLACION DE LAGRANGE 
NO. PUNTOS?6 


INTRODUCIR X(1)?-0.5 
INTRODUCIR Y(1)2?0 
INTRODUCIR X(2)?0.5 
INTRODUCIR Y(2)?0.7 
INTRODUCIR X(3)?2 
INTRODUCIR Y(3)?2 
INTRODUCIR X(4)2?3 
INTRODUCIR Y(4)2?2.6 
INTRODUCIR X(5)?4.5 
INTRODUCIR Y(5)?2.8 
INTRODUCIR X(6)2?6 
INTRODUCIR Y(6)?3.3 


38 PROGRAMAS PARA CIENCIA E INGENIERIA 


NO. PUNTOS EN VENTANA?4 
CALCULAR Y(X),X=?2.3 
Y(2.3)=2.20996 


NUEVA X 1,NUE.VENT2,NUE.EJEC 321 
CALCULAR Y(X), X=?7 
Y(7)=4.4333333 


NUEVA X 1,NUE.VENT2,NUE.EJEC 3?2 
NO. PUNTOS EN VENTANA?5 
CALCULAR Y(X), X=?2.3 
1(2.3)=2,22256 


NUEVA X 1,NUE.VENT 2,NUE.EJEC 321 


CALCULAR Y(X), X=?7 
Y(7)=5.1909091 


NUEVA X 1,NUE.VENT2,NUE.EJEC 321 
CALCULAR Y(X), X=?4.5 


Y(4.5)=2.8 
Figura 9-1 
LEYENDA 
X — MEDIDO 
O - CALCULADO 
y (4 PUNTOS EN VENTANA) 
6 (E - CALCULADO 
(5 PUNTOS EN VENTANA) 
3 O 
10) 


Discusión de resultados 


El ejemplo muestra que la anchura de «ventana» afecta mucho a los va- 
lores extrapolares. El valor interpolado cambió solamente desde 2.20... a 
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2.22..., mientras que el valor extrapolado cambió desde 4.43... a 5.19... El 
usuario ha de decidir qué anchura de «ventana» es más adecuada. El ejem- 
plo muestra también que el polinomio de Lagrange se ajusta exactamente a 
cada punto dado. Cuando Y(4.5) se calcula por el programa, genera 2.8, 
que es el valor exacto originalmente introducido para Y(4.5). 


OBSERVACIONES SOBRE LA PROGRAMACION 


El programa comprueba las entradas correctas en las líneas 80, 170 y 
280. La bifurcación se realiza en tres sentencias simples en las líneas 
560-580. La evaluación de la fórmula de Lagrange tiene lugar en un lazo 
doble en las líneas 430-500. 


LISTADO DEL PROGRAMA 


9 ES 
20 DIM X(20) 
30 DIM Y(20) 
40 PRINT "INTERPOLACION DE LAGRANGE" 
50 PRINT "NO. PUNTOS?"; 
60 INPUT A 
70 PRINT A 
80 IF A<=20 AND A 2>AND A=INT ABS A THEN 
GOTO 110 
90 PRINT "ERROR" 
100 GOTO 50 
110 FOR I=1 TO A 
120 PRINT "INTRODUCIR X("31;")2"; 
130 INPUT X 
140 PRINT X 
150 LET X(I)=X 
160 IF I=1 THEN GOTO 200 
170 IF X(I)>X(I-1) THEN GOTO 200 
180 PRINT "X(1+1)<=X(1)" 
190 GOTO 120 
200 PRINT "INTRODUCIR Y(";31;")2?"; 
210 INPUT Y 
220 PRINT Y 
230, LEE LI=Y 
240 NEXT 1 
250 PRINT "NO. PUNTOS EN VENTANA?"; 
260 INPUT B 
270 PRINT B 
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280 


290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
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IF B<=A AND B>¡ AND B=INT ABS B THEN 


GOTO 310 

PRINT "ERROR" 

GOTO 250 

LET C=INT (B/2) 

LET Es] 

PRINT. “CALCULAR YVOC): X=? 
INPUT W  : 

PRINT Y 


IF W<=X(1) OR I=A THEN GOTO 390 


LETAIsi+1 

GOTO 360 

LET Z=0 

LET D=I-C j 

IF D<l THEN LET D=] 

IF I-C+B>A THEN LET D=A-B+] 
FOR I=D TO D+B-1 

LET E=1 

FOR J=D TO D+B-] 

IF J=1I THEN GOTO 480 

LET E=E*(W-X(3J))/(X(1)-X(3J)) 
NEXT J : 

LEFT Z=Z+E*Y1L1) 

NEXT 1 

PRINT a E a 

PRINT 

PRINT "NUEVA X 1, NUE.VENT.2, 
INPUT F 

PRINT F 

IF F=1 THEN GOTO 310 

IF F=2 THEN GOTO 250 

IF F=3 THEN GOTO 5 

GOTO 530 
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10. Escalas legibles para 
trazados por computadora 


DESCRIPCION DEL PROGRAMA 


Todo buen programa de trazado debe tener la capacidad adecuada para 
calcular automáticamente las escalas para los ejes X, Y o Z. Si los datos a 
trazar se generan por la computadora, el usuario no podrá prever siempre 
las escalas necesarias. Sin embargo, el Timex/Sinclair 1000 puede automa- 
tizar la tarea y, además, hacerlo de forma inteligente, para conseguir 
escalas «legibles». El programa, dado el margen de una variable, produce 
una escala lineal leible con intervarlos uniformes. La salida del programa 
puede utilizarse, entonces, para controlar cualquier rutina de trazado. Una 
escala lineal leible se define aquí como una escala con anchura de intervalo 
igual al producto de una potencia entera de 10 y 1,2 ó 5 y valores de escala 
que son múltiples enteros de la anchura del intervalo. 

La anterior definición permite valores de escala tales como: 


— 0.5, 0.0, 0.5, 1.0.... 
EZ 1,26, 280 
100, 200, 300,..., etc., 


pero prohibiría los ejemplos siguientes: 


Dado un mínimo MIN, un máximo MAX de la matriz a trazar y un nú- 
mero aproximado de intervalos de escala deseados N, el programa calcula 
cuatro parámetros: un nuevo mínimo MINP, un número máximo MAXP, 
una anchura de intervalo DIST y varios intervalos de escala NP (que ajus- 
tarán mejor el trazado y aproximará N. Los parámetros calculados, además 
de producir una escala legible, satisfarán las desigualdades siguientes: 


(MIN — DIST) < MINP<=MIN, 
MAX<=MAXP < (MAX + DIST) 
N/V2.5 < NP < (NV 2.5 + 2). 


Una característica singular del programa es que introduce una puerta 
estrecha =+ J (línea 200) alrededor de los valores MIN y MAX para evitar 
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un margen innecesariamente grande entre MINP y MAXP debido al re- 
dondeo de la computadora. 


INSTRUCCIONES 


Inicializar el programa pulsando la tecla RUN; a continuación seguir 
las indicaciones pata el mínimo y el máximo de su matriz y el número 
aproximado deseado de intervalos de escala. El programa proporciona, en- 
tonces, un nuevo mínimo, un nuevo máximo, las dimensiones del intervalo 
y el número de intervalos. 


EJEMPLOS 


Enunciado del problema 


Calcular los valores de escala para una matriz con MIN, MAX y N, res- 
pectivamente iguales a: 


li = 304, PLE 
2. 345,510, 20 
3. — 1100, — 520, 50 
4. 17.5,17.9, 10 


Ejemplo de ejecución 
PROGRAMA CONSTRUCCION ESCALAS 


INTRODUCIR MINIMO?-3.1 
INTRODUCIR MAXIMO?11.1 

INT. NO. APROX. DE INTERVALOS?5 
NUEVO MIN= -4 

NUEVO MAX= 12 

NUEVO INTERVALO= 2 

NO. INTERVALOS= 8 


INTRODUCIR MINIMO?345 

INTRODUCIR MAXIMO?510 

INT. NO. APROX. DE INTERVALOS?20 
NUEVO MIN= 340 

NUEVO MAX= 510 

NUEVO INTERVALO= 10 

NO. INTERVALOS= 17 
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INTRODUCIR MINIMO?-1100 
INTRODUCIR MAXIMO?-520 
INT. NO. APROX. DE INTERVALOS?50 


NUEVO MIN= - 1100 
NUEVO MAX= =520 
NUEVO INTERVALO= 10 
NO. INTERVALOS= 58 


INTRODUCIR MINIMO?17.9 
INTRODUCIR MAXIMO?17.5 
ERROR, MIN>=MAX 


INTRODUCIR MINIMO?17.5 
INTRODUCIR MAXIMO?17.9 

INT. NO. APROX. DE INTERVALOS?10 
NUEVO MIN= 1249 

NUEVO MAX= 179 

NUEVO INTERVALO=  .05 

NO. INTERVALOS= 8 


Discusión de los resultados 


Los cuatro ejemplos indican que el programa produce escalas «legibles» 
para una amplia gama de entradas. El primer ejemplo daría los valores de 
escala de - 4, - 2, 0, 2,..., 12. El segundo ejemplo proporcionaría 340, 350, 
..., 310, El tercer ejemplo daría como resultado los valores — 1100, — 1090, 
..., 520. El cuarto ejemplo proporcionaría 17.5, 17.55,..., 17.9. El cuarto 
ejemplo muestra también el mensaje de error cuando el mínimo y el máxi- 
mo se introdujeron por error en orden inverso. 


OBSERVACIONES SOBRE LA PROGRAMACION 


La comprobación de errores se produce en la línea 100. La determina- 
ción de las escalas se realiza en las líneas 240-270. La sentencia PRINT en 
la línea 340 realiza también el cálculo para ahorrar código. 


LISTADO DEL PROGRAMA 


20 ¿GL:S 
30 PRINT "PROGRAMA CONSTRUCCION ESCALAS" 
40 PRINT 
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45 

50 

60 

70 

80 

90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
230 
260 
270 
280 
290 
300 
310 
320 
330 
340 


345 
360 
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PRINT "INTRODUCIR MINIMO?"; 
INPUT X 

PRINT X 

PRINT "INTRODUCIR MAXIMO?"; 

INPUT Y 

PRINT Y 

IF Y>X THEN GOTO 130 

PRINT "ERROR, MIN>=MAX" 

GOTO 40 

PRINT "INT. NO. APROX. DE INTERVALOS?"; 
INPUT N 

PRINT ON 

IF N>0 AND N=INT ABS N THEN GOTO 190 
PRINT "ERROR" | 

GOTO 130 

LET D=(Y-X)/N 

LET J=D/1E5 

LET E=INT (LN D/LN 10) 

LET F=D/10**E 

LET V=10 

IF F<SQR 50 THEN LET V=5 

IF F<SQR 10 THEN LET V=2 

IF F<SQR 2 THEN LET V=1 

LET C=V*10**E 

LET G=INT (X/C) 

IF ABS (G+1-X/C)<J THEN LET G=G+1 
LET A=C*G 

LET H=INT (Y/C)+1 

IF ABS (Y/C+1-H)<J THEN LET H=H-1 
LET B=C*H 

PRINT "NUEVO MIN=",A,"NUEVO MAX=",B, 
"NUEVO INTERVALO=", 

PRINT C,"NO. INTERVALOS=",H-G 

GOTO 40 


Originalmente publicado como «Trazados generados por computadora, es- 
calas de calculadora de bolsillo». Reimpresión de Electronics, 15 diciembre 
1981, copyright O McGraw-Hill, Inc 1981. Todos los derechos reservados. 


PROBABILIDAD 
Y ESTADISTICA 


11. Combinaciones, 
permutaciones y factoriales 


DESCRIPCION DEL PROGRAMA 


Este programa evalúa las tres funciones más importantes en estadística 
y en combinatoria. 

Combinaciones. El programa calcula el número de combinaciones posi- 
bles de N objetos diferentes tomados de un conjunto de M objetos. Al 
cambiar la secuencia de los objetos en una combinación no se obtiene una 
nueva combinación. 

La función resultante C(M,N) se denomina también el coeficiente binó- 
mico y puede representarse como sigue: 


Fórmula 11-1 M M! 
_(M+1_ (M1) (Mét_, 
Ñ N a | : l 


La primera expresión, que es la definición del coeficiente binómico 
C(M,N) puede dar lugar a grandes errores al dividir un gran número por 
otro también grande. Si alguna de las factoriales en el numerador o en el 
denominador es mayor que 33!, entonces, no se podría evaluar la primera 
expresión por el Timex/Sinclair 1000. Por el contrario, la segunda expre- 
sión evita el problema del desbordamiento de la capacidad, por lo que será 
la que utilizaremos para evaluar C(M,N). 
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Permutaciones. El programa calcula el número de permutaciones posi- 
bles con M objetos tomados N a N. Al cambiar el orden de los objetos en 
una permutación se obtiene una permutación diferente. 

La función resultante P(M,N) se calcula como sigue: 


Fórmula 11-2 


P (M,N) = M (M — 1)... (M-— N +1) 


Factoriales. Una factorial F(M) es igual al número de formas en las que 
M objetos puede redisponerse. La función se define y calcula como sigue: 


Fórmula 11-3 


F(M) = M (M — 1)....1 


EJEMPLOS 


Enunciados de problemas 


1. ¿En cuántas formas pueden 10 personas disponerse en una fila? Eva- 
luar F(10). 

2. ¿Cuántos grupos de 3 personas pueden seleccionarse de entre un gru- 
po de 20 cuando importa el orden en cada grupo de 3? Evaluar P(20,3). 

3. ¿Cuántos grupos de 3 pueden seleccionarse de entre 20 personas, 
cuando el orden no importa? Evaluar C(20,3). 

4. Lo mismo que en el ejemplo 3 pero seleccionando un grupo de 17 
personas de entre 20. Evaluar C(20,17). 


Ejemplo de ejecución 


FACTORIALES 1,PERMUTACIONES 2 
COMBINACIONES (COEF. BIN.) 321 
FACTORIAL DE?10 
FACT(10)=3628800 


FACTORIALES 1,PERMUTACIONES 2 
COMBINACIONES (COEF. BIN.) 322 
M=?20 

N=? 3 

PERM(20,3)=6840 
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FACTORIALES 1,PERMUTACIONES 2 
COMBINACIONES (COEF. BIN.) 3?3 
M=?20 

N=?3 

COMB(20,3)=1140 


FACTORIALES 1,PERMUTACIONES 2 
COMBINACIONES (COEF. BIN.) 373 
M=?20 

N=?17 

COMB(20,17)=1140 


FACTORIALES 1,PERMUTACIONES 2 
COMBINACIONES (COEF. BIN.) 32 


Discusión de resultados 


El ejemplo 1 produce un número muy grande, los ejemplos 3 y 4 pro- 
porcionan el mismo resultado. La razón para la igualdad es que cuando se 
elige un grupo de 17, el grupo restante será tres, que es lo mismo que en el 
ejemplo 3. 


OBSERVACIONES SOBRE LA PROGRAMACION 


La comprobación de los errores de entrada se realiza en las líneas 60 y 
130. La misma sentencia PRINT en la línea 450 se utiliza para las permu- 
taciones y las combinaciones asignando un nombre adecuado a la variable 
C$. El programa hace uso de la igualdad C(M,N) = C(M,M - N) y seleccio- 
na en la línea 250 el recorrido de cálculo más rápido. 


LISTADO DEL PROGRAMA 


20 LES 

30 PRINT 

35 PRINT "FACTORIALES 1,PERMUTACIONES 2" 
37 PRINT "COMBINACIONES (COEF. BIN.) 32"; 
40 INPUT F 

30 PRENT $ 

60 IF F<>1 AND F<>2 AND F<>3 THEN GOTO 30 
70 GOTO :100*F 

100 PRINT "FACTORIAL DE?"; 

110 INPUT M 

120 PRINT M 
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130 IF M<=0 OR M<>INT ABS M OR M>33 
THEN GOTO 100 

140 LET B=1 

150 FOR I=1 TO M 

160 LET B=B*I 

170 NEXT 1 

180 PRINT "FACT(";M;")=";B 

190 GOTO 30. 

200 LET C$="PERM(" 

210 GOSUB 480 

220 IF M>=N THEN GOTO 250 

230 LET B=0 

240 GOTO 450 

250 LET B=1 ; 

260 FOR I=M-N+1 TO M 

270 LET B=B*I 

280 NEXT 1 

290 GOTO 450 

300 LET C$="COMB(" 

310 GOSUB 480 

320 IF M>=N THEN GOTO 350 

330: LET -B=0 

340 GOTO 450 

350 IF M<>N THEN GOTO 380 

360 LET B=1 

3/70 GOTO 450 

380 LET Y=N8N 

390 IF M<2*N THEN LET Y=M-N 

400 LET K=M+1 

410 LET B=1 

420 FOR I=1 TO Y 

430 LET B=B*(K/I-1) 

440 NEXT 1 

450 PRENT CM. SNS SD 

470 GOTO 30 

480 PRINT "M=2?"; 

490 INPUT M 

495 PRINT M 

500 PRINT "N=2"; 

510 INPUT N 

515 PRINT N 

520 IF M=INT ABS M AND N=INT ABS N AND M*N<>0 
THEN RETURN 

530 PRINT "ERROR" 

540 GOTO 480 


12. La función de error 
ERF y su función ERFC 
complementaria 


DESCRIPCION DEL PROGRAMA 


Los atributos de muchos objetos en la naturaleza, por ejemplo, la altura 
en pulgadas de la población masculina de una cierta ciudad, o los coefi- 
cientes de inteligencia (1Q) de la población, siguen estrictamente la distri- 
bución normal o gausiana (la bien conocida curva de la campana de 
Gauss). Una distribución normal particular se expresa completamente me- 
diante dos parámetros, la media y la desviación cuadrática media. La 
función de error (ERF) es igual a la distribución normal acumulativa y la 
función de error complementaria (ERFC) es igual a 1-ERF según las defi- 
niciones siguientes: 


Fórmula 12-1 
Xx o 
E 15 dt 


Erf (1,0,x) = 
avV2r-—xX 


Las anteriores fórmulas no pueden evaluarse de forma cerrada. El mé- 
todo utilizado en este programa para calcular ERF y ERFC es iterativo y 
lleva a resultados que son correctos dentro de la capacidad de visualización 
completa del Timex/Sinclair 1000 para ambas funciones. Por compa- 
ración, el método normal utilizado para calcular ERF y ERFC se basa en 
la expansión polinómica de un quinto o sexto grado. Y es preciso para so- 
lamente 3 a 4 lugares decimales. La ERFC calculada mediante expansión 
polinómica pierde completamente la exactitud en los extremos, o colas, de 
la distribución. 


INSTRUCCIONES 


Al pulsar la tecla RUN se emiten mensajes de petición de datos para la 
media de distribución, desviación standard y el límite de distribución supe- 
rior en donde deben evaluarse ERF y ERFC. 
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EJEMPLOS 


Enunciados de problemas 


1. Dada una población con un coeficiente de inteligencia (1Q) medio de 
100 y una desviación cuadrática media (desviación estándar) de 15, ¿Qué 
fracción de la población tendrá valores IQ mayores que 135? 

2. El peso medio de un saco de fertilizante es de 45 libras. El proceso de 
fabricación introduce una desviación cuadrática media de 1,5 libras. ¿Qué 
proporción de sacos será más pesada que 49,5 libras? 

3. Probablemente haya oído hablar acerca del límite de sigma de 2,33 
(sigma es la desviación cuadrática media). ¿Qué porcentaje de población 
incluye? 


Ejemplo de ejecución 


PROGRAMA DE FUNCION DE ERROR 
DISTR. MEDIA=?100 

DESVIACION EST.=?15 

LIMITE DISTR. SUPERIOR=?135 
ERF=0.99018467 ERFC=.0098153287 


DISTR. MEDIA=?45 

DESVIACION ST.=?1.5 

LIMITE DISTR. SUPERIOR=?49.5 
ERF=0.9986501 ERFC=.001349898 


DISTR. MEDIA=?100 

DESVIACION ST.=?10 

LIMITE DISTR. SUPERIOR=?123.3 
ERF=0.99009693 ERFC=.0099030756 


Discusión de los resultados 


Solamente el 0,98% de la población tendrá IQ mayores que 135. El se- 
gundo ejemplo indica que el peligro de sacos pesados es bastante pequeño, 
solamente del 0,13%. El tercer ejemplo indica que la desviación estándar 
de 2,33 incluye al 95% de la población. 


OBSERVACIONES SOBRE LA PROGRAMACION 


Dependiendo del valor del argumento, las funciones se evalúan como 
fracciones continuas con un número determinado de elementos o se eva- 
lúan de forma iterativa hasta que convergen. 
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LISTADO DEL PROGRAMA 


¿0 ki) 
30 PRINT 
35 PRINT "PROGRAMA DE FUNCION DE ERROR" 
40 PRINT "DISTR. MEDIA=2"; 
50 INPUT B 
60 PRINT B 
70 PRINT "DESVIACION EST.=?*; 
80 INPUT C 
90 PRINT C 
100 IF C>=0 THEN GOTO 130 
110 PRINT "ERROR" 
120 GOTO 70 
130 PRINT "LIMITE DISTR. SUPERIOR=?"; 
140 INPUT D 
150 PRINT D 
160 LET X=(D-B)/C 
170 LET X=X/SQR 2 
180 LET Y=1/SQR PIXEXP (-X*X) 
190 IF X<SQR 2 THEN GOTO 280 
200 LET A=14/X 
210 FOR 1=27 TO 1 STEP -1 
220 LET A=1/2/(X+A) 
230 NEXT I 
240 LET Z=Y/(X+A) 
250. -LET .25Z/2 
260 LET N=1-Z 
270 GOTO 390 
280 LET I=1 
290 LET T=2*X*Y 
300 LET N=T 
JTAALET PZA RETA L2 LAT) 
320 LET M=N+T 
330 IF M=N THEN GOTO 370 
340 LET N=M 
350 LET I=1+1 
360 GOTO 310 
370 LET N=N/2+.5 
380 LET Z=1-N 
390 PRINT "ERF="3N,*“ERFC="32Z 
400 PRINT 
410 GOTO 40 


13. Distribuciones 
binomia e hipergeométrica 


DESCRIPCION DEL PROGRAMA 


Estas dos distribuciones se producen en muchos problemas que se refie- 
ren a la teoría del muestreo, al control de calidad y a los juegos. La distri- 
bución binomia se aplica al caso de un universo estadístico grande, mien- 
tras que la distribución hipergeométrica se aplica a casos de un universo 
pequeño, en donde las muestras se toman sin ninguna sustitución. Por 
ejemplo, la probabilidad de encontrar X elementos defectuosos en una 
muestra a partir de un universo grande con una proporción conocida de 
elementos defectuosos, seguiría la distribución binomia. Análogamente, la 
probabilidad de encontrar exactamente dos bolas rojas en una muestra de 
cinco bolas tomadas todas a la vez a partir de universo de diez bolas, tres de 
las cuales son rojas, seguiría la distribución hipergeométrica. Ambas distri- 
buciones son muy útiles para estimar acontecimientos de la vida real. Las 
siguientes fórmulas se aplican para las dos distribuciones. 


A. Distribución binomia 


Dado un universo grande con una proporción de P elementos especiales 
(defectuosos, marcados, etc.) y un tamaño de muestra de N, la probabilidad 
de encontrar X elementos especiales en esa muestra es como sigue: 


Fórmula 13-1 


Prob (N,P,X) = (x) y =p 


B. Distribución hipergeométrica. 
Dado un universo (población estadística) de pequeña magnitud de ta- 
maño M con K elementos defectuosos, la probabilidad de encontrar X 
elementos defectuosos en una muestra de tamaño N es como sigue: 


Fórmula 13-2 


Prob (M,N,K,X) = 
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INSTRUCCIONES 


El primer mensaje orientativo después de pulsar la tecla RUN propor- 
ciona la elección de la distribución binomia (1) o hipergeométrica (2). Si se 
selecciona la distribución binomia, sigue una llamada orientativa para el 
tamaño de la muestra N y para la proporción de elementos especiales P en 
el universo. En este y siguientes ejemplos, cuando dicha llamada orientati- 
va exija más de un valor de entrada, ha de introducirse por el teclado cada 
valor por separado seguido por «Enter» (introducir). La siguiente llamada 
corresponde al número de elementos especiales X en la muestra, para lo 
que debe calcularse la probabilidad. Después de visualizar PROB(X) y la 
probabilidad acumulativa, la siguiente llamada orientativa da una elección 
entre una nueva X o una nueva ejecución. 

Análogamente, para la distribución hipergeométrica, la secuencia de 
llamadas orientativas es la misma, con la excepción de que se requiere la 
introducción de tres parámetros M, K y N (tamaño del universo, número 
de elementos especiales y el tamaño de la muestra) para describir la distri- 
bución. 


EJEMPLOS 


Enunciados de problemas 


1. Distribución binomia. Se sabe que el 4% de la marca de un fabrican- 
te es de elementos defectuosos. ¿Cuál es la probabilidad de encontrar, a lo 
sumo, un elemento defectuoso en una muestra de cinco? 

2. Distribución hipergeométrica. En una caja de 20 chocolates, 5 son 
oscuros y 15 son claros. ¿Cuál es la probabilidad de encontrar uno o más 
chocolates oscuros en un puñado de tres? 


Ejecuciones de ejemplos 


BINOMIA 1,HIPERGEOMETRICA 2?] 


INTRODUCIR N,P?5 . 04 

X=? 1 

PROB(1)=0.16986931 

SUMA DE P(X)= 0.16986931 
NUEVA X 1,NUEVA EJEC. 2?1 

Xx=?0 

PROB(0)=0.8153727 

SUMA DE P(X)= O 98524201 


NUEVA X 1, NUEVA EJEC. 2?2 


BINOMIA 1,HIPERGEOMETRICA 2?2 
INTRODUCIR M,K,N 20 
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5 3 
X=? 1 
PROB(20,5,3,1)= 0.46052632 


SUMA DE (X)=0.46052632 


NUEVA X 1,NUEVA DISTR. 2?1 

X=?2 

PROB(20,5,3,2)= 0.13157895 
SUMA DE (X)=0.59210526 

NUEVA X 1,NUEVA DISTR. 2?1 

X=73 

PROB(20,5,3,3)= .0087719298 
SUMA DE(X)=0.60087719 

NUEVA X 1,NUEVA DISTR. 2?1 


Xx=?0 
PROB(20,5,3,0)=0.39912281 
SUMA DE(X)=1 

NUEVA X 1,NUEVA DISTR. 2? 


Discusión de los resultados 


En el primer ejemplo, la probabilidad de encontrar un elemento defec- 
tuoso en una muestra es 0.169..., y de encontrar uno sin defecto es 0.815... 
La probabilidad de encontrar 0 ó 1 elementos defectuosos en una muestra 
de 5 es, pues, P(0) + P(1) = 0,985... 

En el segundo ejemplo, la probabilidad de encontrar uno o más choco- 
lates oscuros en una muestra de tamaño tres es la suma de probabilidades 
de encontrar uno, dos o tres y es igual a 0.6. El mismo resultado puede ob- 
tenerse más fácilmente calculando la probabilidad de ningún chocolate 
oscuro y restándola de 1.0, como se muestra en la última ejecución del 
ejemplo. 


OBSERVACIONES SOBRE LA PROGRAMACION 


Una comprobación amplia de los valores de entrada se realiza en las lí- 
neas 90, 150, 190, 280, 360 y 400. Los coeficientes binomios se calculan en 
la subrutina entre las líneas 520 y 620 de la misma forma que en el progra- 
ma 11. Obsérvese el empleo de sentencias PRINT en las líneas 230 y 460, 
que incluyen la reimpresión de parámetros de entrada para una más fácil 
interpretación de los resultados. 
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LISTADO DEL PROGRAMA 


20 
40 
41 
45 
50 
60 
70 
80 
90 
100 
110 
120 
130 
140 
150 


160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
390 
360 


365 
370 


380 
330 


CLS 

PRINT 

PRINT "BINOMIA 1,HIPERGEOMETRICA 2?"; 
LET Z=0 

INPUT Y 

PRINT <u 

IF W=1 THEN GOTO 100 

IF W=2 THEN GOTO 290 

GOTO 40 

PRINT "INTRODUCIR N,P2"; 

INPUT N 

PRINT ON, 

INPUT P 

PRINT P | 

IF N<>INT N OR N<=0 OR P<=0 OR P>=1 
THEN GOTO 100 

PRINT "X=2"; 

INPUT X 

PRINT X 

IF X<>INT X OR X<0 OR X>N THEN GOTO 160 
GOSUB 520 | 

LET T=B*P**x*(]-P)**(N-X) 

LET Z=Z+T 

PRINT "PROB(";X3")=";T,"SUM DE P(X)=",Z 
PRINT "NUEVA X 1,NUEVA EJEC.?2"; 
INPUT $4 

PRINT Qu 

IF W=1 THEN GOTO 160 

GOTO 40 

PRINT "INTRODUCIR M,K,N", 

INPUT E 

PRINT E, 

INPUT G 

PRINT G, 

INPUT F 

PRINT F 

IF E<>INT E OR G<>INT G OR F<>INT F 
THEN GOTO 290 

IF E<=0 OR F<=0 OR G<=0 OR F>E OR G>E 
THEN GOTO 290 

PRINT "X=2"; 

INPUT D 

PRINT D 
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400 IF D<>INT D OR D<0 OR D>G THEN GOTO 370 
410 LET N=G 

420 LET X=D 

422 GOSUB 520 

424 LET J=B 

426 LET N=E-G 

428 LET X=F-D 

430 GOSUB 520 

432 LET J=J*B 

434 LET N=E 

436 LET X=F 

438 GOSUB 520 

440 LET J=J/B 

450 LET Z=Z+J : 

460 PRENT "PROB("SES". Tb. 
465 .PRINT FS" 5D70)=*% 0, "SUM. DECA]S"SZ 
470 PRINT "NUEVA X 1,NUEVA DISTR. 2?"; 
480 INPUT W 

490 PRINT W 

500 IF W=1 THEN GOTO 370 

510 GOTO 40 

520 LET Y=X 
530 IF X>(N-X) THEN LET Y=N-X 
540 IF Y<>0 THEN GOTO 570 

550 LET B=1 

560 RETURN 

570 LET B=1 

580 LET L=N+] 

590 FOR I=1 TO Y 

500. .LET. BsB*(L1=1) 

610 NEXT I 

620 RETURN 


14. Desviantes 
pseudo-aleatorias 
(gaussianas) normales 
y uniformes 


DESCRIPCION DEL PROGRAMA 


Este programa produce una secuencia de desviantes aleatorias (gaussia- 
nas) normalmente distribuidas o uniformes. Dichas secuencias suelen utili- 
zarse para la realización de modelos con el método de Monte Carlo, en 
donde no sea factible un cálculo directo. 

Un ejemplo de una aplicación para la que no es factible un cálculo 
exacto es el de la satisfacción de los abonados telefónicos. La satisfacción 
de un grupo de abonados telefónicos depende de la distribución probabilís- 
tica de la sonoridad de quienes hablan, de la distribución del ruido y de la 
distribución de pérdidas de transmisión en una conexión telefónica. Aun- 
que cada uso de estas distribuciones puede estimarse como gaussiana con 
una desviación estándar y media conocida a partir de encuestas, la determi- 
nación de la distribución combinada es difícil desde el punto de vista del 
cálculo. Sin embargo, el efecto combinado de todos estos factores puede si- 
mularse fácilmente generado, de forma repetida, desviantes pseudo-alea- 
torias a partir de cada una de estas CRU IONES y multiplicando las 
probabilidades individuales. 

Este programa actuaría mamita como una subrutina en una apli- 
cación proporcionada por el usuario. Las fórmulas utilizadas para generar 
las desviantes pseudo-aleatorias normales y uniformes son las siguientes: 


A. Distribución uniforme 


Fórmula 14-1 
X, = Fract (997 X ,, + 7) 


en donde RND es la función pseudo-aleatoria del Timex/Sinclair 1000 
BASIC y X e Y son los límites inferior y superior de la variable aleatoria. 
La anterior fórmula genera los valores de distribución uniforme entre los 
límites inferior y superior X e Y. 
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B. Distribución normal (gaussiana) 


Fórmula 14-2 
X, = (— 2 In U)” sin (27 U,) 


en donde U, y U, son dos desviantes pseudo-aleatorias uniformemente dis- 
tribuidas, y secuencialmente generadas, entre O y 1. Las desviaciones pseu- 
do-aleatatorias normalmente distribuidas resultantes tienen una media de O 
y una desviación estándar de 1. Los valores se transforman, entonces, por 
el programa en una distribución con una media arbitraria M y una desvia- 
ción estándar S multiplicando cada desviante aleatoria con S y añadiendo 
luego M. 


INSTRUCCIONES 


Pulsar la tecla RUN para iniciar el programa y luego, seguir las llama- 
das orientativas para desviantes uniformes (1) o gaussianas (2). La siguiente 
llamada orientativa es para el elemento de arranque que inicializa la fun- 
ción aleatoria (RND) del Timex/Sinclair 1000. El valor del elemento de 
arranque debe estar comprendido entre O y 65535. Al introducir O se ini- 
ciará siempre una secuencia aleatoria diferente; de no ser así, la secuencia 
se determinará previamente. A continuación, seguir las llamadas orientati- 
vas para los límites inferior y superior para la distribución uniforme o para 
la desviación estándar o media de la distribución gaussiana, respectiva- 
mente. 


EJEMPLOS 


Enunciados de problemas 


1. Generar varias desviantes uniformemente distribuidas entre 10 y 20. 

2. Generar desviantes gaussianas con una media de 100 y una desvia- 
ción estándar de 10. Una posible interpretación de los valores resultantes 
de las desviantes gaussianas podría ser una lectura aleatoria de los cocien- 
tes intelectuales (IQ). 


Ejemplos de ejecución 


UNIFORME 1, GAUSSIANA 2?1 
ELEMENTO ARRANQUE=?123 

LIMITE INFERIOR, SUPERIOR?10,20 
DESVIANTE ALEATORIA=11.418915 


MAS VALORES 1, NUEVA EJEC. 221 
DESVIANTE ALEATORIA=16.428375 
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MAS VALORES 1, NUEVA EJEC. 2?1 
DESVIANTE ALEATORIA=12.132111 


MAS VALORES 1, NUEVA -EJEC.* 271 
DESVIANTE ALEATORIA=19.917297 


MAS VALORES 1, NUEVA EJEC. 2?1 
DESVIANTE ALEATORIA=13.797302 


MAS VALORES 1, NUEVA EJEC. 2?] 
DESVIANTE ALEATORIA=14.804688 


MAS VALORES 1, NUEVA EJEC. 2?2 


UNIFORME 1, GAUSSIANA 2?2 
ELEMENTO ARRANQUE=?555 
MEDIA=?100 

DEV. “Edd. =7Y 10 

DESVIANTE ALEATORIA=90.646835 


MAS VALORES 1, NUEVA EJEC. 2?] 
DESVIANTE ALEATORIA=114.96507 


MAS VALORES 1, NUEVA EJEC. 2?1 
DESVIANTE ALEATORIA=108.74587 


MAS VALORES 1, NUEVA EJEC. 2?1 
DESVIANTE ALEATORIA=116.17717 


MAS VALORES 1, NUEVA EJEC. 2?1 
DESVIANTE ALEATORIA=100.46652 


MAS VALORES 1, NUEVA EJEC. 2?1 
DESVIANTE ALEATORIA=114.59181 


Discusión de los resultados 


Las desviantes uniformes están comprendidas entre 10 y 20 y las des- 
viantes gaussianas siguen la distribución normal. 


OBSERVACIONES SOBRE LA PROGRAMACION 


La comprobación de errores de todas las entradas se realiza en las líneas 
130, 220 y 280. El programa muestra el empleo de las funciones RAND y 
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RND de BASIC. Si el valor del elemento de arranque, o argumento, de 
RAND se selecciona como 0, una nueva secuencia aleatoria se iniciará 
cada vez que se ejecute el programa. De cualquier otro modo, los valores 
pseudo-aleatorios pueden repetirse, si se elige el mismo elemento de arran- 
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que, en ejecuciones sucesivas. 


LISTADO DEL PROGRAMA 


20 
25 
30 


CELS 

PRINT 

PRINT "UNIFORME 1, GAUSSIANA 2?"; 
INPUT B 

PRINT B ] 

PRINT "ELEMENTO ARRANQUE=?; 

INPUT S 

PRINT S 

RAND (5) 

IF B=1 THEN GOTO 80 

IF B=2 THEN GOTO 230 

GOTO 30 

PRINT "LIMITE INFERIOR,. SUPERIOR?"; 
INPUT X 

PRINT X; 

INPUT Y 

PRINTS EV 

IF Y>X THEN GOTO 160 

PRINT "MIN>=MAX" 

GOTO 80 

PRINT "DESVIANTE ALEATORIA=";RND*(Y-X)+X 
PRINT 

PRINT "MAS VALORES 1, NUEVA EJEC. 2?”; 
INPUT Z 

PRINT Z 

IF Z=1 THEN GOTO 160 

IF Z=2 THEN GOTO 25 

GOTO 170 

PRINT. "MEDIAS ?2"= 

INPUT C 

PRINT C 

PRINT "DEV. EST.=?"; 

INPUT D 

PRINT D 

IF D>=0 THEN GOTO 310 

PRENDA DUE Ve. Est 0) 


300 
310 
315 
320 
20 
330 
340 
350 
360 
370 
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GOTO 255 

PRINT "DESVIANTE ALEATORIA="; 

PRINT SQR (-2*LN RND)*SIN (2*PI*RND)*D+C 
PRINT 

PRINT "MAS VALORES 1, NUEVA EJEC. 22?"; 
INPUT Z 

PRINT Z 

IF Z=1 THEN GOTO 310 

IF Z=2 THEN GOTO 25 

GOTO 320 


MATEMATICAS 


15. Coeficientes Fourier 
de funciones periódicas 


DESCRIPCION DEL PROGRAMA 


Una función periódica continua y = f(x) con un período T definido en- 
tre - T/2 y +T/2 puede representarse como la suma de términos de las 
funciones seno y coseno, según las siguientes fórmulas clásicas atribuidas 
al matemático francés Fourier. 


Fórmula 15-1 


- A S EL T ,_ As : ul 
f(x) = 3 + 2 (A cos2n q x + B, sin 2n 5 Xx) = o QnA x+0) 
T/2 
A, = 4 100 cos 72 ax C, =4/ A?+ B? 
-T/2 Ñ A 
B, = 2 Fo nr 9 = Arctan(- Pl, 
E T An 


-T/2 


Este programa evalúa las anteriores expresiones realizando una integra- 
ción numérica a intervalos de muestreo especificados. Si se conoce, por 
anticipado, que la función es par, f(x) = f(- x), o que es impar, f(x) =- f(- x), 
el programa puede acelerarse porque sólo han de evaluarse los coeficientes 
A o B y la integración puede proseguir en solamente una mitad del período 
T. El programa se aprovecha de esta información. 

Aunque el programa solamente visualiza los coeficientes de Fourier su- 
cesivos, estos coeficientes podrían almacenarse fácilmente en el momento 
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del cálculo en una variable dimensionada (que permita la capacidad de la 
memoria) para su uso futuro. 


INSTRUCCIONES 


Antes de ejecutar el programa, ha de escribirse una subrutina que defi- 
na la función a evaluar comenzando en la línea 1000 y finalizando con 
RETURN (ver listado del programa y ejemplo). Al pulsar la tecla RUN se 
inicializa el programa y se emite una llamada orientativa para el tipo de 
función: (1) para una función par, (2) para una función impar y (3) si no se 
conoce el tipo de función. Las llamadas orientativas posteriores requieren 
la magnitud del período T y el número de intervalos de muestreo (dos o 
más). La selección de un número grande de intervalos proporciona una 
mayor exactitud pero requiere más tiempo de cálculo. Para una función 
par o impar, el número requerido de intervalos se aplica a solamente una 
mitad del período T, con lo que se hacen más pequeños los intervalos y 
mejor la exactitud del cálculo. Cuando se determina un coeficiente de Fou- 
rier, la computadora indicará su valor de forma intermitente. Los sucesivos 
coeficientes se calculan y visualizan hasta que la pantalla esté llena (error 
(5) o se pulse la tecla BREAK. 


EJEMPLO 


Enunciado de problema 


La función especificada en las líneas 1000-1020 se muestra en la figura 
15-1. Su periodo está comprendido entre — 3.5 y + 3.5. El programa calcula 
los coeficientes Fourier sucesivos AO, A1, A2... para esta función par de x. 
La curva aproximada se muestra a trazos en la misma figura. El cálculo ha 
de realizarse con 10 y 100 intervalos. 


Figura 15-1 
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Ejemplo de ejecución 


FUNCION PAR 1, IMPAR 2, 
NINGUNA DE ELLAS 3?1] 


PERIODO T=?7 


NO. INTERVALOS=?10 


A0=4.025 
Al=3.0177238 
A2=1.1368327 
A3=0.1117739 
A4=.070788988 
A5=0.14849242 
A6=.058580645 
A7=.01478557 
A8=.027038986 
A9=.011989939 


FUNCION PAR 1, IMPAR 2, 
NINGUNA DE ELLAS 321 


PERIODO T=?7 


NO. INTERVALOS=?100 


A0=4.00015 
A1=3.0345888 
A2=1.1793983 
A3=0.10390619 
A4=.058466219 
A5=0.18842451 
A6=.08414665 
A7=.0001493158 
A8=.047250562 
A9=.057932911 


A10=0 A10=.0094089601 
A11=-.011989941 A11=.0077683475 
A12=-.027038986 A12=.03243459 
A13=-.014785568 A13=.017811636 
A14=-.058580644 A14=.00014726251 
A15=-0.14849242 A15=.013311727 
A16=-.070788986 A16=.018110545 
A17=-0.1117739 


Discusión de los resultados 


Los coeficientes de Fourier se verificaron de forma independiente con 
un algoritmo de alta precisión en una computadora diferente; se encontró 
que con 100 intervalos entre cuatro y cinco lugares decimales son correctos 
en el Timex/Sinclair 1000. Obsérvese el fenómeno de Gibbs en esquinas 
agudas, X =- 3.5, 0, +3.5; la expansión de Fourier no converge en esos 
puntos a su valor verdadero. Por ejemplo, a X =0 la serie parece converger 
hacia 6,85 en lugar de a 7,0 como estaba previsto. En otros puntos, la con- 
vergencia es excelente. 


OBSERVACIONES SOBRE LA PROGRAMACION 
Las comprobaciones de los errores de entrada se realizan en las líneas 


60, 100 y 140. El programa realiza la integración trapezoidal de la función 
comenzando en la línea 1000. 
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LISTADO DEL PROGRAMA 


30 PRINT "FUNCION PAR 1,IMPAR 2,NINGUNA DE 
ELLAS ARAS 

40 INPUT A 

50 PRINTAA 

60 IF A<] OR A>3 OR A<>INT A THEN GOTO 30 

70 PRINT *PERITODO T=72*> 

80 INPUT B 

90 PRINT B 

00 IF B<=0 THEN GOTO 70 

O PRINT "NO. INTERVALOS=?"; 

O INPUT C 

O PRINT C 

O IF C<2 OR C<>INT C THEN GOTO 110 

O IF A=3 THEN GOTO 500 

O LET D=B/2/C 

O IF A=2 THEN GOTO 340 

180 LET J=-1 

190 LET J=J+1 

195 PAUSE 20 

200 LET X=0 

210 GOSUB 1000 

220 LET E=Y 

230 LET F=0 

240 LET G=2*J*PI/B 

CD OR: JT OE 

260 LET X=D*I 

270 GOSUB 1000 

280 LET F=F+(E+Y)/2*D*C0S .(6*(X7D42J ) 

290 LET E=Y 

300 NEXT 1 

310 LET H=4*F/B 

3200 PRINT. “Aud 3510 

330 GOTO 190 

340 LET J=0 

350 LET J=J+1 

390 PAUSE>:20 

360 LET X=0 

3/0 GOSUB 1000 

380 LET E=Y 

390 LET K=0 

400 LET G=2*J*PI/B 

410 FOR I=1 TO C 

420 LET X=D*I 
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430 GOSUB 1000 

440 LET K=K+(E+Y)/2*D*SIN (G*(X-D/2)) 
450 LET E=Y 

460 NEXT 1 

470 LET L=4*K/B 

480. PRINT "B"5J3*="3L 

490 GOTO 350. 

500 LET D=B/C 

510. LET.JS=1 

520 LET J=J+1 

525 PAUSE. 20 

530 LET X=-B/2 

540 GOSUB 1000 

550 LET E=Y 

560 LET K=0 

570 LET F=0 

580 LET G=2*J*PI/B 

590 FOR I=1 TO C 

600 LET X=D*I-B/2 

610 GOSUB 1000 

620 LET M=G*(X-D/2) 

630 LET F=F+(E+Y)/2*D*C0S M 
640 IF J=0 THEN GOTO 660 
650 LET K=K+(E+Y)/2*D*SIN M 


660 LET E=Y 

670 NEXT 1 

680 LET H=2*F/B 

690 PRINT "A";J;"=";H 


700 IF J=0 THEN GOTO 520 
710 LET L=2*K/B 

720 PRENT "B"30;"="3L 

730 LET N=SQR (H*H+L*L) 

740 LET P=ATN (-L/H) 

750: BRINT “ON 

760 PRINT "THETA";J;"=";P 
770 GOTO 520 

1000 LET Y=7-3.5*ABS X 

1010 IF ABS X>2 THEN LET Y=0 
1020 RETURN 


16. Soluciones de ecuaciones 
diferenciales por el método 
de Runge-Kutta 


DESCRIPCION DEL PROGRAMA 


Los métodos de Runge-Kutta para soluciones numéricas de ecuaciones 
diferenciales ordinarias están constituidos por un grupo de fórmulas iterati- 
vas que se distinguen por algunas características comunes importantes. 

1. La derivada y” puede darse como una función de x y de y. No ha de 
proporcionarse como una función explícita de x solamente. 

2. El método es de iniciación automática. El conocimiento de la fun- 
ción sólo se requiere en un solo punto para calcular el punto sucesivo en la 
iteración. 

3. El intervalo de integración puede cambiarse cuando prosigue la ite- 
ración. : 

4. Las fórmulas suelen converger rápidamente. 

Sucesivamente, fórmulas de Runge-Kutta de más alto orden e interva- 
los de integración más pequeños dan lugar a soluciones más precisas. El 
programa evalúa las fórmulas de Runge-Kutta de cuarto orden para resol- 
ver ecuaciones diferenciales de primer orden del tipo: 


y' = f(x,y) 


En las fórmulas se supone que se conoce f(x,y). La condición inicial de 
X0 y YO, el intervalo de integración H a lo largo del eje X y el valor final 
de X se seleccionan por el usuario. El programa evalúa las fórmulas si- 
guientes: 


Fórmula 16-1 
K, = Hf(X,, Y;) 
K, = Hf(X, + H/2, Y, + K//2) 


K, = Hf(X +H/2, Y +K./2) 


67 


68 PROGRAMAS PARA CIENCIA E INGENIERIA 


K, = Hf(X, + H/2, Y, +K.y/2) 


Y = Y, + (K,+2K,+2K,+K,) 


INSTRUCCIONES 


Antes de ejecutar el programa, escribir la función y = f(x,y) como una 
subrutina de BASIC que comienza en la línea 1000 y añadirla al programa 
como se muestra en el listado del programa y en el ejemplo. Utilicese la le- 
tra R para la derivada. Al pulsar la tecla RUN se inicia el programa y se 
produce la llamada orientativa para la magnitud del intervalo H, los valo- 
res iniciales de X0 y YO y para el valor final de X en el que debe encontrar- 
se la solución. El intervalo H debe ser positivo, aunque X puede ser más 
pequeña o más grande que X0. No hay otras limitaciones para H. El pro- 
grama hace que parpadeen los valores x sucesivos a medida que se evalúan 
y se para cuando se termina el cálculo visualizado Y(X). 


EJEMPLO 


Enunciado del problema 


Dada la ecuación diferencial y' =x + y, con la condición inicial X =0 y 
Y =0. Determinar Y(- 1). Suponer los intervalos de integración de 0,15 y 
de 0,05. 


Ejemplo de ejecución 


RUNGE KUTTA DE CUARTO ORDEN 
INTERVALO DE INTEGRACION H?0.15 
INTRODUCIR X0,YO,XLAST?O0,1 
0:18:0.:3,0.45:50,6,0.75,0.95 1 
Y(1)=0.71827251 


INTERVALO DE INTEGRACION H?.05 
INTRODUCIR XO,YO,XLAST?0,0,1 
-05,0.-5,0.139,0:2,0.290:,0:3.0,39 
AO A550. 500.357 00507092: U0 8 
1550:8:0.85,0.:9,0,.9% 515 
Y(1)=0.71828169 


¿0 
0. 
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Discusión de resultados 


Se trata de una ecuación diferencial para la que puede encontrarse una 
solución exacta. La solución es: 


Y(1) = e - 2 = 0.718281... 


La solución Runge-Kutta es 0,718272... para H=0,15 y 0,718281 para 
el intervalo de integración de 0,05. El más corto intervalo de integración 
proporciona evidentemente una solución más exacta. 


OBSERVACIONES SOBRE LA PROGRAMACION 


La comprobación de los errores de entrada sobre la magnitud del inter- 
valo se realiza en la línea 60. Obsérvese la visualización del valor final de 
Y(X) en la línea 245. En la sentencia se utiliza XLAST como su argumento 
para una mayor claridad. Como se indicó en la descripción del programa, 
la rutina que comienza en la línea 1000 describe la ecuación diferencial 
particular objeto de evaluación. 


LISTADO DEL PROGRAMA 


20 PRINT "RUNGE KUTTA DE CUARTO ORDEN" 
30 PRINT "INTERVALO DE INTEGRACION H?"; 
40 INPUT A 

50 PRINT A 

60 IF A<=0 THEN GOTO 30 

70 PRINT "INTRODUCIR XO,YO,XLAST?"; 

80 INPUT B 


90 PRINT B; 

100 INPUT C 

110 PRINT *,*50; 

120 INPUT D 

130 PRIKT 2D 

135 IF B=D THEN GOTO 70 
140 LET E=INT ABS (D-B)/A 
150 LET F=B 

160 LET G=SGN (D-B)*A 

170 LET H=C 


180 FOR I=1 TO E 

190 GOSUB 270 

200 NEXT 1 

210 IF F=D THEN GOTO 240 
220 LET G=D-F 

230 GOSUB 270 
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240 
245 
¿ou 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
1000 
1010 


PRINT 

PRINT O o 
PRINT 

GOTO 30 

LET X=F 

LET Y=H 
GODUB 1000 
LET K=G*R 
LET X=F+G6/2 
LET Y=H+K/2 
GOSUB 1000 
LET L=G*R 
LET Y=H+L/2 
GOSUB 1000 
LET M=G*R 
LET. XA7F+G 
LET Y=H+M 
GOSUB 1000 
LET N=G*R 
LET H=H+K/6+L/3+M/3+N/6 
LET F=F+6G 
PRINTS 
PAUSE 20 
RETURN 

LET R=X+Y 
RETURN 


17. Soluciones numéricas 
de ecuaciones 
trascendentes 


DESCRIPCION DEL PROGRAMA 


Este programa resuelve ecuaciones numéricas que, por lo general, no 
pueden resolverse fácilmente de forma explícita. Por ejemplo, resolverá 
para X: 


F(X) = LOGX - X*+ 10=0. 


El programa solicita una estimación inicial de X0O en las proximidades 
de la solución. Dicha estimación puede hacerse habitualmente a partir de 
un trazado de la función. A continuación, solicita un incremento DELTA 
en el que se evaluará la función (un DELTA grande convergerá más rápi- 
damente, pero puede omitir la solución verdadera). Ahora, el programa 
solicita EPSILON, un pequeño número positivo. Cuando el valor absoluto 
de la función es menor que EPSILON, se interrumpe la iteración y se im- 
prime la solución. Si existe más de una solución, el programa encontrará la 
solución más próxima a la estimación inicial. 

El programa utiliza un método de iteración de Newton modificado para 
encontrar el cruce por cero de la función, que es la solución requerida. El 
programa comienza por «trazan» una línea recta a través de la función eva- 
luada en X0 + DELTA/2 y X0 - DELTA/?. El cruce de esta línea con el 
eje X da la siguiente estimación de XO. Si la estimación de X0 interrumpe 
el cambio mientras el valor absoluto de F(X) siga siendo más grande que 
EPSILON, el intervalo DELTA se dividirá por la mitad por el programa y 
continuará la iteración. Es evidente que este método solamente actuará 
para funciones sin irregularidades entre la primera estimación de XO y el 
valor verdadero de X0. La función F(X) ha de ser también de valor único y 
su primera derivada no puede ser O en esta zona; de cualquier otro modo, 
la línea horizontal recta no cruzará el eje X. Para indicar cómo la iteración 
sigue adelante (icon la esperanza de que haya convergencia!), el programa 
hace parpadear las aproximaciones sucesivas de X0. En tanto que estas úl- 
timas sean convergentes, debe llegar a la solución aproximada en un breve 
período de tiempo. Entonces, el programa se interrumpe y visualiza su me- 
jor estimación de la solución correspondiente a la ecuación. 
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INSTRUCCIONES 


Como en el programa anterior, la ecuación a resolver se añade al pro- 
grama antes de su ejecución. La subrutina que describe la ecuación debe 
comenzar en la línea 1000 y ha de finalizar con RETURN, tal como se in- 
dica en el listado y ejemplo. Al pulsar la tecla RUN se inicia el programa y 
genera tres mensajes orientativos, para la conjetura inicial de X0, DELTA 
y EPSILON. A continuación, la iteración continúa hasta que se verifique 
ABS(F(X)) < EPSILON. 


EJEMPLO 


Enunciado del problema 


Trataremos de resolver la ecuación antes citada 


F(X) = LOGX -— X?* + 10 =0. 


Utilicemos la conjetura inicial de X0=2.0, DELTA =0.01 y EPSI- 
LON = 1E - 6 y luego, se repetirá con una DELTA más grande de 0.1. 


Ejemplo de ejecución 


SOLUCION ECUACION TRANSCENDENTE 
INTROD. CONJ. PARA XO,DELTA X,EPS.? 


ds 0d DE =0 

APROXIMACIONES SUCESIVAS: 
3.6656754 LO 13I0I3 
3.24213/4 3.2420352 
0 LA 2zU se SL ALUOZ 
3.2420352 3.2420382 
3.2420382 3.242039 
3.242039 3.2420392 


X0=3.2420392 


SOLUCION ECUACION TRANSCENDENTE 
INTROD. CONJ. PARA XO,DELTA X,EPS.? 


2:03 LE, TESÓ 

APROXIMACIONES SUCESIVAS: 
3,.6650052 3.266929 

3.2417394 3.2416374 
3.2416374 3.2416374 


SOLUCIONES NUMERICAS DE ECUACIONES TRASCENDENTES 


.2419388 
.2419388 
.2420141 
.2420329 
2420311 
.2420388 
.2420391 


0D 0 0 UY UY UY) Uy 


X0=3.2420391 


Discusión del problema 


0 0 0 UU UM UU 


.2419388 
.2420141 
.2420141 
.2420329 
.2420377 
.2420388 
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Ambas ejecuciones, después de parpadear los valores convergentes de 
X, conducen al resultado de X = 3.242... 


OBSERVACIONES SOBRE LA PROGRAMACION 


El flujo del programa finaliza con STOP en la línea 260 para no interfe- 
rir con la subrutina que describe la ecuación. Para cada ejecución, se debe 


pulsar la tecla RUN 


LISTADO DEL PROGRAMA 


10 PRINT "SOLUCION ECUACION TRANSCENDENTE" 
PARA XO DELTA Xx EPS¿2?” 


20 PRINT "INT. 
30 INPUT A 
35 LET E=A 
40 PRINT A; 
50 INPUT B 


60 PRINT.“,"GB>s 


70 INPUT C 
SO PRENTO"S 736 


85 PRINT "APROXIMACIONES SUCESIVAS:" 


90 LET X=A-B/2 
100 GOSUB 1000 
110 LET D=Y 
120 LET X=A+B/2 
130 GOSUB 1000 


CONJ. 


140 LET Z=B*D/(D-Y)+A-B/2 


150 LET X=Z 
160 PRINT Z, 
170 PAUSE. 30 
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180 GOSUB 1000 

190 IF ABS Y<=ABS C THEN GOTO 240 
200 IF Y=E THEN LET B=B/2 
210 LET A=Z 

220 LET:E=sY 

230 GOTO 90 

240 PRINT 

245 PRINT 

250 PRINT "X0="3Z 

260 STOP 

1000 LET Y=LN X/LN 10-X*X+10 
1010 RETURN 


18. Aritmética 
de números complejos 


DESCRIPCION DEL PROGRAMA 


Con este programa, su computadora Timex/Sinclair 1000 se hace una 
calculadora aritmética de números complejos sofisticada que trabaja en la 
notación polaca inversa (RPN). Para quienes estén familiarizados con las 
calculadoras Hewlett Packard, que utilizan RPN, ésta será anticuada. Si no 
se está familiarizado con RPN se necesitarán algunas explicaciones. 

Consideremos una pila de registros, con cuatro en profundidad, deno- 
minados X, Y, Z y T en su computadora. La parte inferior de la pila, el 
registro X y el situado inmediatamente encima, el registro Y, siempre serán 
objeto de visualización. Cada registro contiene las partes real e imaginaria 
de un número complejo. Las operaciones usuarias (de una sola variable), 
como son la elevación al cuadrado, la extracción de raíz cuadrada, el cálcu- 
lo de recíprocos, etc., se registran solamente en el registro X y no afectan a 
la pila. Las operaciones con dos variables, tales como las de suma, resta, 
multiplicación y división, trabajan en los registros X e Y. El resultado de 
una operación de dos variables aparece siempre en el registro X y descen- 
derá la pila; esto es, el contenido del registro Z se desplaza al registro Y y el 
registro T se copia en Z. El registro T conserva su valor anterior. Los nú- 
meros se introducen en la pila mediante la orden ENTER (1). Al introducir 
un número complejo se le desplaza al registro X sustituyendo el contenido 
anterior de ese registro. Se proveen varias Operaciones de manipulación de 
la pila: X<—-> Y intercambia el contenido de los registros X e Y, RAISE 
STACK destruye el contenido anterior del registro T, desplaza el conteni- 
do anterior del registro Z a T, de Y a Z y duplica X en Y. El programa 
proporciona también un registro de memoria M. Se dispone de dos órdenes 
para copiar el registro X en M (SAVE) y para copiar el registro M en X 
(RECALL). 


INSTRUCCIONES 


RUN — Comenzar el cálculo, emitir llamadas orientativas para las opera- 
ciones l a 12 
1 —-—ENTER, introduce el número complejo en el registro X y la pila 
se mantiene invariable 
2 — —RAISE STACK, copia el contenido de Z en T, de Y en Z y de X 
en Y; destruye el contenido anterior de T 
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3 —X<-—->Y, intercambia los registros X e Y 
4 — "SAVE, copia X en M 
5 —HRECALL, copia M en X 
6 — Suma X e Y, pone el resultado en X y hace bajar la pila 
7 — Resta X de Y, pone el resultado en X y hace bajar la pila 
8 — Multiplica X e Y, pone el resultado en X y hace bajar la pila 
9 — Divide Y por X, pone el resultado en X y hace bajar la pila 
10 — Calcula el valor recíproco de X y lo pone en X 
11 ——Calcula el valor de la raíz cuadrada de X y lo pone en X 
12 — Calcula el cuadrado de X y lo pone en X 


EJEMPLO 


Enunciado del problema 


Evaluar la expresión siguiente: 


Fórmula 18-1 


Q + 35) (5 — 43) 
5 + 4j 


Ejemplo de ejecución 


OPERACION 1-12?1 
INTRODUCIR R,1?2,3 


REAL IMAGINARIA 
Y: 0 0 
ie 3 
OPERACION 1-12?3 

REAL IMAGINARIA 
eZ 3 
X: 0 0 


OPERACION 1-12?1 
INTRODUCIR R,1?5,-4 


REAL IMAGINARIA 
Va 2 3 
ob -4 
OPERACION 1-12?8 

REAL IMAGINARIA 
Y> “UY 0 


ES LL / 


ARITMETICA DE NUMEROS COMPLEJOS 77 


OPERACION 1-12?3 


REAL IMAGINARIA 
YE RZ 7 
Xx: 0 0 


OPERACION 1-12?1 
INTRODUCIR R,I?5,4 


REAL IMAGINARIA 
ME ZA 7 
X: 5 4 
OPERACION 1-12?9 

REAL IMAGINARIA 
0 0 
A: 33058537 -1.2926829 
OPERACION 1-12?11 

REAL IMAGINARIA 
vo. Q 0 
X: 1.8670036 -0.34619187 
OPERACION 1-12?12 

REAL IMAGINARIA 
Ys 9 0 
XK 34300037 -1.2926829 


Discusión de resultados 


Se introduce primero el número complejo 2+3j. A continuación se 
desplaza desde el registro X al Y antes de que se introduzca el segundo nú- 
mero complejo 5-4j. Luego, se multiplian X e Y y el resultado inmediato es 
22 +7j. A continuación el producto se desplaza al registro Y y se introduce 
el último número complejo 5 +4j. Después de esta operación, se realiza 
una división y una raíz cuadrada. El resultado es, entonces: 1.86... 
- 0.346... j. En el paso final, el resultado se eleva al cuadrado para demos- 
trar que la raíz cuadrada y la elevación al cuadrado de números complejos 
son Operaciones inversas. Obsérvese que existe también una segunda raiz 
que es la negativa de la primera, a saber: 

- 1.86...+ 0.346...j. 


OBSERVACIONES SOBRE LA PROGRAMACION 


El programa llena la capacidad de RAM de 2K del Timex/Sinclair 
1000. Por consiguiente, sólo los registros X e Y se visualizan en lugar de 
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la pila completa. La función raíz cuadrada (11) visualiza una sola raíz. La 
otra raíz es la negativa de la primera. La elevación y descenso de la pila se 
realiza en los lazos FOR NEXT en las líneas 100-110 y 1000-1015. La úni- 
ca comprobación de los errores de entrada se realiza en la línea 27. La 
bifurcación múltiple se realiza en la línea 40 multiplicando el código de 
operación por 50. 


LISTADO DEL PROGRAMA 


10 DIM X(8) 
15 PRINT 

18 PRINT  “"OPERACTON: T=T28%3 

ZO. INPUT+ FE É 

25 PRINT F 

27. TE FS=INT FOR. FS] ORIE=TZ2 TREN -G0TO 18 
40 GOTO 50*F 

50 PRINT "INTRODUCIR R,I1?"; 

55 INPUT X(1) 

60 PRINT X(1); 
65 INPUT X(2) 
67 PRINT *,* 
( 
( 


73 PRINT " REAL","IMAGINARIA","Y: 
EE E e 

74 LET D=RSR+JI*J 

23: GOTO 15 

00. .FOR, F=8: TO 3 STEP -1 
LET X(1)=X(1-2) 


180 GOTO 70 
200 LET K=R 
205 LET M=J 
210 GOTO 70 
250 LETS 


300 LET X(1)=R+X(3) 
305 LET X(2)=J+X(4) 


310 
350 
30) 
360 
400 
405 
410 
450 
455 
460 
500 
505 
510 
550 
591 
90% 
553 
554 
09 
50U 
600 
605 
610 
1000 
1010 
1015 
1020 


ARITMETICA DE NUMEROS COMPLEJOS 


GOTO 1000 

LET ACTIEXLI)=R 

LET X(2)=X(4)-J 

GOTO 1000 

LET X(1)=R*X(3)-J*X(4) 

LET X(2)=R*X(4)+J*X(3) 
GOTO 1000 

LET X(1)=(R*X(3)+J*X(4))/D 
LET X(2)=(R*X(4)-J*X(3))/D 
GOTO 1000 

LET X(1)=R/D 

LET X(2)=-3J/D 

GOTO 1000 

IF R<>0 OR J<>0 THEN GOTO 554 
LET X(1)=0 

LET X(2)=0 

GOTO 70 


LET X(1)=SQR ((R+SQR (R*R+J*J))/2) 
LET X(2)=J/2/X(1) 

GOTO 70 

LET X(1)=R*R-J*J 

LET X(2)=2*R*J 

GOTO 70 

FOR I=3 TO 6 

LET X(1)=X(1+2) 

NEXT 1 

GOTO 70 
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19. Paso de fracciones 
decimales 
a fracciones binarias 


DESCRIPCION DEL PROGRAMA 


La parte del mundo civilizado que emplea el sistema británico de uni- 
dades de media todavía prefiere”las fracciones binarias tales como 5/16 6 
7/32 a los números tan fáciles de recordar como 0.3125 ó 0.21875. Para 
atender a ambos sistemas, este programa convertirá una fracción decimal 
en su equivalente en fracción binaria. Definiremos una «fracción binaria» 
como una fracción cuyo denominador es una potencia entera de 2. Todo lo 
que el usuario del programa tiene que hacer es indicar la más pequeña frac- 
ción aceptable con un denominador potencia binaria de 2, v.g. 256 por 
1/256, 64 por 1/64, etc. 

Si la conversión puede hacerse con exactitud, el programa dividirá el 
numerador por el denominador. Por ejemplo, 0.375 se convertirá a 3/8, 
aunque la más pequeña fracción requerida por el usuario puede ser 1/256. 
Si la conversión no puede realizarse con exactitud, el más pequeño deno- 
minador binario se encontrará dentro de las limitaciones dadas y se 
indicará el error porcentual. 


INSTRUCCIONES 


Al pulsar la tecla RUN se emiten dos llamadas de mensaje orientativo: 
INTRODUCIR FRACCION DECIMAL? (el número a convertir en la 
fracción binaria y DENOMINADOR MAX. (POT. DE 2) (el número a 
indicar el mayor denominador admisible). Entonces, el programa visualiza 
la fracción binaria y el error, si lo hubiere. 


EJEMPLO 


Enunciado del problema 


Convertir 2:5, 0.12717 y 0.1875 en fracciones binarias. No aceptar nin- 
guna fracción más pequeña que 1/16 para el primer ejemplo y 1/256 para 
los ejemplos segundo y tercero. 
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Ejemplo de ejecución 


INTRODUCIR FRACCION DECIMAL?2.5 
DENOMINADOR MAX (POT. DE 2)2?16 
2.955/2 


INTRODUCIR FRACCION DECIMAL?0.12717 
DENOMINADOR MAX. (POT. DE 2)?256 
0.12717 APROXIMA 33/256 

CON ERROR DE -1.3652984 POR CIENTO 


INTRODUCIR FRACCION DECIMAL?0.1875 
DENOMINADOR MAX. (POT. DE 2)?256 
0.1875=3/16 


Discusión de resultados 


Los ejemplos primero y tercero conducen a fracciones binarias exactas. 
El segundo ejemplo puede aproximarse dentro de las limitaciones de 1/256 
para dentro de 1.35...%. 


OBSERVACIONES SOBRE LA PROGRAMACION 


La única comprobación de errores con respecto al valor de entrada (de- 
nominador más grande) se realiza en las líneas 70-90. Para aplicar el 
redondeo del Timex/Sinclair 1000 se utiliza un factor «ficticio» de 0.0001 
en las líneas 90 y 120. Para encontrar el más pequeño denominador, el 
programa hace que se divida el numerador y el denominador por 2 (líneas 
140-150), en tanto que ambos sean enteros. 


LISTADO DEL PROGRAMA 


5 PRINT 
10 PRINT "INTRODUCIR FRACCION DEC.?"; 
20 INPUT A 
30 PRINT A 
40 PRINT "DENOMINADOR MAX. (POT. DE 2)2?"; 
50 INPUT B 
60 PRINT B 
70 IF B<2 THEN GOTO 40 
80 LET C=LN B/LN 2 
90 IF ABS (C-INT C)>.0001 THEN GOTO 40 
100 LET C=INT € 
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130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
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LET D=A*B 

IF (ABS D-INT ABS (D+.0001)>=.0001) 
THEN GOTO 170 

IF D/2<>INT (D/2) THEN GOTO 230 

LET D=D/2 

LET B=B/2 

GOTO 130 

LET D=INT (D+.5) 

LET E=A-D/B 

LET F=100*E/A 

PRINT A;" APROXIMA ";D;"/";B 

PRINT "CON ERROR DE ";F;" POR CIENTO" 
GOTO 5 

PRINT Ate SE 

GOTO 5 


INVESTIGACION 
OPERATIVA 


20. Solución de ecuaciones 
de las colas 


DESCRIPCION DEL PROGRAMA 


Este programa resuelve las ecuaciones de las colas M/M/S que descri- 
ben muchas situaciones de la vida real. 

Estas ecuaciones suponen una distribución exponencial de las tasas de 
llegada de clientes L, una distribución exponencial de los tiempos de servir 
M, S sirvientes idénticos y un orden de servicio FIFO (primero en entrar, 
primero en salir). La distribución exponencial implica independencia de 
eventos; los clientes buscarán servicio con independencia de la longitud de 
la cola y los sirvientes trabajarán a un régimen constante haciendo caso 
omiso de la carga y de la longitud de la cola. Dichos supuestos suelen pro- 
porcionar estimaciones seguras (pesimistas) de los tiempos de espera y de 
otros parámetros de colas. Los casos de M/M/S dan lugar a ecuaciones re- 
lativamente simples (en términos de la teoría de las colas) y pueden 
resolverse dentro de las limitaciones de memoria del Timex/Sinclair 1000. 
El programa puede utilizarse para estimar problemas tales como cuántos 
sillones de barbero hay que proporcionar para un establecimiento situado 
en el centro de la ciudad, cuántas operadoras de teléfonos deben trabajar 
entre las 2 y las 3 de la tarde o cuántas impresoras deben conectarse a un 
procesador central. Por lo general, la aplicación del programa ayudará a 
compensar el coste de proporcionar servicio con la incomodidad de largas 
colas de espera. Las ecuaciones de las colas utilizadas en el programa se 
dan a continuación. 


Limitaciones del programa: S < 34, U < =1.0, todas las entradas posi- 
tivas, todos los resultados intermedios deben ser inferiores a 1.0E38. 


Entradas: 
Tasa de llegada media r 
Tasa de servicio media m 
Número de sirvientes s 
Números de elementos en sistemas n 
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Salidas: 
Utilización de las instalaciones U = r/(ms) 
Prob. de no encontrar ningún elemento en sistema P(0) 


P(0) =((1/m)' /sI1-U)+ E (my 507 


Prob. de encontrar n elementos en el sistema P(n) 


P(n) =P(0) x (r/m)" x (1/n!) da n<s 
P(n) =P(0) x (r/my x (1A(s!s"=)) P ESP 


Prob. de encontrar ocupados a todos los sirvientes B 


B= Ss P(n) = (1/m)y xP(OM(s! x (1-U)) 


Tiempo de espera medio en cola Tw = B/(sm x (1 —- U)) 

Tiempo de respuesta medio del sistema Tr = Tw+ 1/m 

Dev. estánd. de Tr=Tw/Bx V  (BQ-B)+s? x (1 — U)?) 

Prob. de Tw > T=BxExp -smT x (1 - U)) 

Prob. de Tr > T =Exp (-mT) x (1 +(1 —- Exp - msTk)) x B/sk) para 
k< > 0 en donde k=1-U-I/s 

Prob. de Tr > T=Exp -mT)x(1 +BmT) para k=0 

Número medio de elementos en cola Q =r x Tw 

Número medio de elementos en sistema N =r x Tr 


El programa proporciona todas las llamadas de mensajes orientativos y 
comprueba los datos de entrada para ver si los valores son correctos. 


INSTRUCCIONES DEL PROGRAMA 


Al pulsar la tecla RUN emitirá llamadas orientativas para introducir la 
tasa de llegada L, la tasa de servicio M y el número de servicio S. Una en- 
trada incorrecta hará que se repitan las llamadas orientativas. A continua- 
ción, el programa calcula y visualiza los parámetros anteriormente citados 
del sistema de las colas, a saber B, U, P(0), Tw, Tr, desviación standard de 
Tr, Q y N. Entonces, el programa emitirá un mensaje orientativo: «IN- 
TRODUCIR COD 1, 2 OR 39». Para calcular P(n), introducir 1,P(Tw, 
Tr)- 2 o para continuar con una nueva ejecución — 3. Al seleccionar P(n) se 
generará un mensaje orientativo para el número de elementos en el sistema 
n y calculará y visualizará P(n). Al elegir P(Tw,Tr) se emitirá un mensaje 
para el tiempo T y calculará y visualizará P(Tw > T) y P(Tr >T). A con- 
tinuación, se repiten los mensajes orientativos. 


EJEMPLO 


Enunciado del problema 


Una zona de trabajo de computadora tiene cinco terminales. Por tér- 
mino medio, 16 clientes por hora llegan a intervalos aleatorios con cada 
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trabajo teniendo una desviación media de 15 minutos (cuatro por hora). 
Calcular los parámetros de las colas y encontrar también la reducción en el 
tiempo de espera si se instala un sexto terminal. 

Ejemplo de ejecución 


INTRODUCIR L,M,S?16,4,5 


P(TODO OCUPADO)= 0.55411255 
UTIL.= 0.8 

P(0)= 012987013 
ESPERA MED.= 0.13852814 
RESP.MED.= 0.38852814 
DESV.ESTAND.(TR)= 0133392003 
MED. EN SIST.= 6.2164502 
MED. EN COLA= 2.2164502 


INTRODUCIR CODIGO 1,2 O 321 
NO. ELEM. EN SISTEMA?] 
P(N=1)=.051948052 


INTRODUCIR CODIGO 1,2 O 3?] 
NO. ELEM. EN SISTEMA?6 
P(N=6)=.08865809 


INTRODUCIR CODIGO 1,2 O 3?1 
NO. ELEM. EN SISTEMA?10 
P(N=10)=.03631432 


INTRODUCIR CODIGO 1,2 O 322 
TIEMPO=?0.2 
P(TW>0.2)=0.24897882 
P(TR>0.2)=0.64851202 


INTRODUCIR CODIGO 1,2 O 3?2 
TIEMPO=?10 
P(TW>10)=2.3540663E-18 
P(TR>TF0)=9,84110.1:2E-17 


INTRODUCIR L,M,S?16,4,6 
P(TODO OCUPADO)= 0.28476085 


UTIL.= 0.66666667 
P(0)= .016685206 
ESPERA MED.= .035595106 
RESP.MED. = 0. ¿855951] 


DEV.ST.(TR)= 0.26482403 
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MED. EN SIST.= 4.5695217 
MED. EN COLA= 0.56952169 


INTRODUCIR COD. 1,2 O 3?1 
NO. ELEMENTOS EN SIST.?1 
P(N=1)=.066740823 


INTRODUCIR CODIGO 1,2 O 32?2 
TIEMPO=?0.2 
P(TW>0.2)=.057492223 
P(TR>0.2)=0.51978804 


Discusión de resultados 


La primera parte del ejemplo (cinco terminales) muestra una utiliza- 
ción bastante alta de las instalaciones, 0.8. El tiempo de espera medio en 
cola es 0.318 o algo más de 8 minutos. La adición de un sexto terminal 
producirá una mejora considerable. El tiempo de espera disminuirá a 0.035 
horas o unos dos minutos. Todas las demás probabilidades de espera resul- 
tarán también afectadas de forma correspondiente. 


OBSERVACIONES SOBRE LA PROGRAMACION 


Una comprobación de los valores de entrada se realiza en la línea 80. 
Obsérvese el empleo del factor ficticio de redondeo en la línea 60. Una sen- 
tencia PRINT bastante complicada se muestra en la línea 900. A medida 
que se introduzca el programa se observará el empleo de la función EXP 
(desplazar ENTER para introducir el modo FUNCION, y luego pulsar la 
tecla X). 


LISTADO DEL PROGRAMA 


9 LS 

10 PRINT "INTRODUCIR L,M,S?"; 

20 INPUT L 

30 PRENTOES 

40 INPUT M 

50: PRINT "SM; 

60 INPUT S 

TOIPRENE AS 

80 IF L<=0 OR M<=0 OR S<=0 OR S<>INT S OR S>33 
THEN GOTO 10 

85 LET P=0 


90 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
ecu 
230 
240 
250 
260 
270 


£ 15 
280 


290 


350 
355 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
910 
520 
530 
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LET X=L/M 

LET U=X/S 

LET V=1-U 

IF V<=0 THEN GOTO 10 
LET Y=S 

GOSUB 1000 

LET T=Z 

FOR I=0 TO S-1 

LET Y=I 

GOSUB 1000 

LET. PEPEX**I/Z 

NEXT 1 

LET Y=S 

GOSUB 1000 

LET P=ILL(PFAFXAS/T/MV) 
LET BEE SEP TLV 

LET W=B/S/M/V 


LET R=W+1/M 

PRINT "P(TODO OCUPADO)=",B,"UTIL=",U, 
"Pp 1013, 

PRINT P,"ESPERA MED.=",W,"RESP.MED.=",R 


PRINT "DESV.ESTAND.(TR)=",SQR(2*B-2*B+(S*V) 
**2)*W/B , 

PRINT “MED. CEN SIST.=".L*R, "MED: 
EN COLA=",L*W 

PRINT 

PRINT "INTRODUCIR CODIGO 1,2 O 32"; 
INPUT K 

PRINT K 

IF K<>1 THEN GOTO 510 

PRINT "NO. ELEM. EN SISTEMA?"; 
INPUT N 

PRINT N 

LETD=P*X*EN 

IF N<S THEN GOTO 460 
LET“ESDLTASRNIN=S) 

GOTO 490 

LET Y=N 

GOSUB 1000 

LET E=D/Z 

PRINT "P(N="3N3")=";E 

GOTO 350 

IF K<>2 THEN GOTO 5 

PRINT "TIEMPO=2"; 

INPUT H 
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540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
900 
910 
1000 
1010 
1020 
1030 
1040 


PRINT H 

LET F=B*EXP -(S*M*V*H) 

LET A$S="4"" 

GOSUB 900 

LET A$="R" 

LET G=V-1/S 

LET D=EXP -(M*H) 

IF ABS G<=.00001 THEN GOTO 640 
LEFIF=SDF(IRBLSELISEXP <(M*S*G*H) 1/6) 
GOTO 650 

LET F=D*(1+B*M*H) 

GOSUB 900 

GOTO 350 

PRINT PATAS PSA 
RETURN 

LET Z=1 

FOR J=1 TO Y 

LET Z=Z*J 

NEXT y 

RETURN 


Originalmente publicado como «La calculadora de bolsillo aborda los pro- 
blemas clásicos de las colas». Reimpresión de Electronics, 22 septiembre 
1981, copyright O McGraw-Hill, Inc. 1981. Todos los derechos reservados. 


21. Parámetros 
de fiabilidad 


de sistemas complejos 


DESCRIPCION DEL PROGRAMA 


La fiabilidad de un sistema compuesto (v. g. central telefónica constitui- 
da por enlaces interurbanos de entrada y de salida, dispositivos de conmu- 
tación y fuentes de alimentación o un centro de cálculo constituido por 
unidades CPU, impresoras, unidades de disco y fuentes de alimentación) 
puede expresarse en función del Tiempo medio antes del fallo (MTBF) y 
del Tiempo medio hasta la reparación (MTTR). Un alto MTBF y un bajo 
MTTR puede conseguirse utilizando componentes fiables que no se dete- 
rioren con frecuencia y sean fáciles de reparar o proporcionando redun- 
dancia. Por ejemplo, dos unidades centrales de proceso (CPU) pueden 
compartir la carga; si una falla, la otra ocupará su lugar. 

El programa determinará la fiabilidad total de un sistema dados los pa- 
rámetros MTBF y MTTR de componentes individuales. Se supone que los 
fallos de componentes individuales son independientes entre sí. La curva 
de fiabilidad sigue, entonces, la parte central de la curva de «pila de baño» 
típica, en donde el alto número de fallos al principio es causado por pro- 
blemas iniciales y el alto número de fallos hacia el final de la vida útil es 
causado por la degradación por el uso del componente. También se supone 
que los componentes individuales se agruparán en subsistemas de uno de 
los tres tipos citados a continuación y que varios subsistemas están conec- 
tados en serie para formar un sistema total. 

A. Subsistemas en serie. El fallo de un solo componente incapacitará el 
sistema completo. 

B. Subsistemas en paralelo. Al menos un componente ha de mantener- 
se activo para que el sistema pueda funcionar. 

C. Subsistema K de n en paralelo. En donde K componentes de entre n 
han de estar en buen estado (activos) para que el sistema funcione. 

El programa calcula MTBF y MTTR para cada subsistema y para el 
sistema total. También determina el factor de disponibilidad y la probabi- 
lidad de que x componentes fallen durante el tiempo t, en donde pueden 
seleccionarse x y t. 

Se utilizan las fórmulas siguientes: 
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A—Disponibilidad 

U-No disponibilidad 

N—n.? de elementos en subsistema 

k—n.* de elementos requeridos para la explotación del sistema 
ñ—n.* de fallos durante el tiempo t (valor medio) 

P(k,t)-Prob. de k fallos durante el tiempo t 

¡i—<lemento, s-subsistema, subíndice 


Fórmulas generales: 


Ñ MTBF o 
- MTBF + MTTR á n = MTBF 
—ñ n k 

Subsistemas en serie: 

N a | 1—-A 
A, = TA, MTBF, = 1/2 (1/MTBF) MTTR, = 2er MTBF, 

1 1 s 
Subsistemas en paralelo: 

N A 

A, = 1 — T (1 — A) MTBF, = Í A MTTR, 


N 
MTTR, = 1/2 (1/MTTR) 
1 


«K de N» en paralelo 


MTTR, = MTTR (N-—K+1) 


MTBE 1... (N—K)! (K=D! 
MTBF, = MTBF, IS 


MTTR, 
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INSTRUCCIONES 


Pulsar la tecla RUN para iniciar el programa y luego, seguir el mensaje 
orientativo introduciendo un número 0 y 6. 


0 — Comenzar una nueva ejecución 

| — Introducir MTBF y MTTR para un subsistema en serie 

2 — Introducir MTBF y MTTR para un subsistema en paralelo redundante 

3 — Introducir MTBF y MTTR para un subsistema k de n en paralelo re- 
dundante 

4 — Borrar subsistema si se cometió un error durante la introdución 

5 — Incluir un subsistema en un sistema total 

6 — Calcular la probabilidad de x fallos durante el tiempo t. 


EJEMPLO 


Enunciado del problema 


Considérese una instalación de central telefónica, tal como se muestra 
en la figura 21-1, que está constituida por una fuente de alimentación y 
otra de reserva, tres sistemas de conmutación, dos de los cuales se requie- 
ren para el funcionamiento, un grupo de enlaces interurbanos de entrada y 
otro de salida conectados a antenas de microondas, todo con parámetros 
MTBF y MTTR conocidos. Calcular los parámetros de fiabilidad total in- 
cluyendo la probabilidad de fallo en 12 y en 24 horas. 


Figura 21-1 
FUENTE DE SISTEMAS ENTRADA/SALIDA 
ALIMENTACION CONMUTACION (EN SERIE) 
EN PARALELO (2 de 3) 


ALIMENTACION 
COMERCIAL 


MTBF= 480 
MTTR= 2 


ENLACES 
INTERURBANOSP 
ENTRADA 


ENLACES 
INTERURBANOS 
SALIDA 


MTBF= 100 
MTTR= 10 


MTBF= 5000 MTBF= 4000 
MTTR=0.5 MTTR=0.5 


ALIMENTACION 
RESERVA 


MTBF= 100 
MTTR=2 
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Ejemplo de ejecución 


OPERACION O - 6722 
EN PARALELO 
INTRODUCIR MTBF,MTTR?480,2 


A= A. IITBADBE 
MTBF = 480 
MTTR= ' 2 


OPERACION O - 6?2 
EN PARALELO 
INTRODUCIR MTBF,MTTR?100,10 


A= 0.99962278 
MTBF= 4416.666/ 
MTTR= 1.6666667 
OPERACION O - 625 

ASYS= 0.-.99962278 
MTBFSYS= 4416.666/7 
MTTRSYS= 1.6666677 


OPERACION O - 623 

K-0F-N 

INTRODUCIR MTBF,MTTR?100,2 
INTRODUCIR K,N?2,3 


A= 0.99880144 
MTBF= (sde BE RO A ES 
MTTR= 1 
OPERACION O - 625 

ASYS= 0.99842467 
MTBFSYS= 701.0582 
MTTRSYS= 1.1061378 
OPERACION O - 62] 

EN SERIE 

INTRODUCIR MTBF,MTTR?5000,0.5 
A= 0.99990001 
MTBF= 5000 

MTTR= 0.50000025 


OPERACION O - 621 

EN SERIE 

INTRODUCIR MTBF,MTTR?4000,0.5 
A= 0.99977504 
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= 0.99977504 
MTBF= ¿222.222 
MTTR= 0.50002885 
OPERACION O - 625 
ASYS= 0.99820007 
MTBFSYS= IZ SUL 
MTTRSYS= 0.96097057 


OPERACION O - 676 
TIEMPO,NO.FALLOS?24,1 
P(24,1)=.043050892 


OPERACION O - 626 
TIEMPO.NO.FALLOS?12,1 
P(12,1)=4022015632 


OPERACION O - 626 
TIEMPO,NO.FALLOS?24,2 
P(24,2)=.00096937613 


Discusión de resultados 


Al examinar el subsistema, vemos primero que el subsistema de fuentes 
de alimentación en paralelo tiene un MTBF de 4416 horas, el subsistema 
de dispositivos de conmutación de 2 de entre 3 tiene un MTBF de 833 ho- 
ras y el subsistema de enlaces interurbanos de entrada/salida en serie tiene 
un MTBF de 2222 horas. El sistema total tiene un MTBF de 532 horas. 
Para mejorar este factor se requeriría evidentemente una mejora en la fiabi- 
lidad de los dispositivos de conmutación. La probabilidad de un solo fallo 
en 24 horas es 4,3%, de un solo fallo en 12 horas es 2,2%. La probabilidad 
de 2 fallos en 24 horas es solamente de 0.097%. 


OBSERVACIONES SOBRE LA PROGRAMACION 
El programa llena casi la totalidad de la memoria de 2K del Timex/Sin- 


clair 1000: Por este motivo, los mensajes orientativos tuvieron que reducir- 
se al mínimo. 
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LISTADO DEL PROGRAMA 


15 GOSUB 800 

20 LET M=1 

40 LET N=0 

50 PRINT 

60 PRINT "OPERACION O - 62"; 
70 INPUT Z 

80 PRINT Z 

90 GOTO 100*Z+10 

10 PRINT "EN SERIE" 

O GOSUB 1000 

0. LET CsCEILA 

O LET D=D*A/(A+B) 

O EET ST 

O LET H=F*(1-D)/D 

O LET E=D 

O PRINT "A=",E,"MTBF=",F,"MTTR=",H 
O GOTO 50 

210 PRINT "EN PARALELO" 
215 GOSUB 1000 

220 LET C=C+1/B 

225 LET D=D*B/(A+B) 

¿30 LET ¿E=1=D 

240 LET H=1/C 

290.LET E=E*H/D 

260 GOTO 180 

3101 ARINT “¡"“K=-0F-N*" 

315 GO0SUB 1000 

320 PRINT "INTRODUCIR K,N?"; 
325 INPUT K 

330 PRINT K; 

390 LNPUT- R 

380 PRINTS 3R 

345 LET L=R-K 

350 GOSUB 900 

355 LET F=J 

360 LET L=K-] 

365 GOSUB 900 

370 LET F=F*J 

3/5 LET L=R 

380 GOSUB 900 

385 LET F=A*F/J*(A/B)**(R-K) 
390 LET H=B/(R-K+1) 

395 -LET: EsF/(F+H) 
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400 GOTO 180 

410 GOSUB 800 

420 PRINT "SUBSISTEMA BORRADO" 
430 GOTO 50 

510 GOSUB 800 

520 LET M=M*E 

530 LET N=N+1/F 

540 LET P=(1-M)/M/N 

550 LET X=1/N 

560 PRINT "ASYS=";M,"MTBFSYS=",X,' MTTRSYS=",P 
570 GOTO 50 

610 PRINT "TIEMPO,NO.FALLOS?*; 
612 INPUT T 

6.14 PRINTS 

615 INPUT Q 

620 PRINT. 5350 

625 LET S=T*N 

630 LET L=Q0Q 

635 GOSUB 900 

640 "LETMUSEXP SESPROF 

645 PRINT dd ar eo ED E 
650 GOTO 50 

800 LET C=0 

810 LET D=1 

820 RETURN 

900 LET J=1 

910 FOR I=1 TO L 

920 LEFT J=3*1 

930 NEXT 1 

940 RETURN , 

1000 PRINT "INTRODUCIR MTBF,MTTR2"; 
1010 INPUT A 

+OZ2O0MPRINT A 

1030 INPUT B 

1040 PRINT ",";B 

1050 RETURN 


MISCELANEA 


22. Frecuencia de las notas 
musicales en la escala 
bien templada 


DESCRIPCION DEL PROGRAMA 


Este programa determina la frecuencia de una nota musical dada en la 
escala bien afinada. La escala supone una relación de frecuencias constante 
de la raíz 12-ésima de 2 entre dos notas contiguas. Johann Sebastián Bach 
estaba tan excitado con la idea de la escala bien templada (afinada) que es- 
cribió una serie de composiciones de piano en la clave de cada nota en esa 
escala. Las octavas en el programa están numeradas como enteros negati- 
vos O positivos, o 0, con la nota A («la») en la octava O con la frecuencia de 
440 Hz. La frecuencia de las notas se dobla con cada octava sucesivamente 
más alta. 


INSTRUCCIONES 


Al pulsar la tecla RUN se inicia el programa y se emiten dos mensajes 
orientativos correspondientes a la nota y al número de la octava. La nota 
debe introducirse como C, C+, D, D+, E, F, F+,G, G+,A, A+o Bi; la 
octava como — 2, -—1,0, 1,2... etc. Como el Timex/Sinclair 1000 no tiene 
una tecla de sostenido (+ ), este signo se sustituye por (+). 


EJEMPLO 


Enunciado del problema 


Determinar las frecuencias de la nota A («la») en la octava 0, C+ en la 
tercera octava, D en la octava — 4 y C y B en la primera octava (1). 


(1) Nota del T.: Las equivalencias de la nomenclatura anglosajona de las notas musicales 
es la siguiente: 

C-«do»; D-«re»; E-«mi»; F-«fa»; G-«sol»; A-«la», B-«si»; C+ =«do» sostenido; D+= «re» 
sostenido; F+= «fa» sostenido; G+ = «sol» sostenido; A+= «la» sostenido. 
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Ejemplo de ejecución 


INTRODUCIR NOTA E.G. C,C+,D..?A 
INTRODUCIR OCTAVA?0 
A,OCTAVA 0=440 HZ 


INTRODUCIR NOTA E.G, C,C+,D..?C+ 
INTRODUCIR OCTAVA?3 
C+,OCTAVA 3=2217.4611 HZ 


INTRODUCIR NOTA E.G. C,C+,D..?D 
INTRODUCIR OCTAVA? -4 
D,OCTAVA -4=18.354048 HZ 


INTRODUCIR. NOTA E.Gi CsbtiDe 2 K 
INTRODUCIR OCTAVA?3 
NOTA ERRONEA 


INTRODUCIR NOTA E.G. Ci, C+,D..?2U 
INTRODUCIR OCTAVA?1 
C,OCTAVA 1=523.25113 HZ 


INTRODUCIR NOTA E.G. C,C+,D..?B 
INTRODUCIR OCTAVA?1 
B,OCTAVA 1=987.7666 HZ 


Discusión de resultados 


El primer ejemplo, como se esperaba, proporciona 440 Hz, la cuarta 
entrada muestra una introdución errónea (K) que proporciona un retorno 
de error y la repetición del mensaje orientativo. 


OBSERVACIONES SOBRE LA PROGRAMACION 


En este caso, para cambiar, se trata de un programa corto. Las líneas 
100 y 135 realizan la comprobación de errores de las entradas. Si la desig- 
nación de las notas no se reconociera como legítima, entonces, la variable I 
mantiene su valor asignado de 100 y se visualiza la observación (comenta- 
rio) «NOTA ERRONEA». El cálculo de la frecuencia se realiza en la línea 
170. 
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C,C+,D. 


LISTADO DEL PROGRAMA 
10-UES 
50 PRINT 
52 PRINT "INTRODUCIR NOTA E.G. 
60 INPUT B$ 
70 PRINT B$ : 
80 PRINT "INTRODUCIR OCTAVA?*; 
90 INPUT B 
100 IF B<>INT B THEN GOTO 80 
FO PRINT B 
120 LET I=100 
121 IF B$="C" THEN LET I=1 
122 IF B$="C+" THEN LET I1=2 
123 IF B$="D" THEN LET I=3 
124 IF B$="D+" THEN LET I=4 
125 IF B$="E" THEN LET I=5 
126 IF B$="F" THEN LET I=6 
122 IF BH="EF*" THEW-LET ==? 
128 IF B$="G" THEN LET I= 
129 IF B$="G+" THEN LET T=9 
130 IF B$="A" THEN LET I=10 
131 1 BS="4+*" THEN LE "111 
132 IF B$="B" THEN LET I=12 
135 IF 1<>100 THEN GOTO 170 
150 PRINT "NOTA ERRONEA" 
160 GOTO 50 
170 LET F=2**((1-10)/12+B)*440 
180 -PRINT 'B$;", OCTAVA "3B3"="3F;3" HZ" 
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GOTO 50 


23. Juego de estrategia 
de inteligencia artificial 


DESCRIPCION DEL PROGRAMA 


Este programa, que es una versión del popular juego de CARA o CRUZ 
es bastante singular por cuanto que aprende jugando con el oponente hu- 
mano. De este modo, demuestra un fuerte rasgo de inteligencia. 

El juego consiste en que el jugador elige CARA o CRUZ. Si la compu- 
tadora tiene una respuesta correcta al hacer la misma elección, entonces 
ganará la computadora. Si la elección es diferente, ganará el jugador. Si us- 
ted, el jugador, juega al azar, la computadora también lo hará en sus 
intentos de acierto y el número de veces en ganar o perder será casi el mis- 
mo. Sin embargo, si sigue alguna de varias estrategias simples; v.g. jugar 
siempre CARA o jugar siempre CRUZ hasta que pierda, y entonces cam- 
biar a CARA y viceversa, o jugar lo mismo (CARA o CRUZ) si pierde, 
pero de forma alternada si gana, la computadora adivinará su estrategia 
después de unas pocas jugadas y le derrotará todas las veces desde entonces, 
hasta que cambie su estrategia. Tan pronto como cambie su estrategia y 
pierda la computadora, esta última comenzará a jugar al azar hasta que 
«adivine» su nueva estrategia y le vuelva a ganar. 

El programa comprueba las dos jugadas anteriores en el juego y consul- 
ta su memoria para el mismo tipo de situación anteriormente jugada. Si se 
encuentra una situación anterior, la computadora hará un intento de acier- 
to basándose en la experiencia anterior. Si entonces, el jugador actúa según 
lo previsto, ganará la computadora. Si la situación no se produjo con ante- 
rioridad, la computadora jugará al azar y almacenará el resultado de su 
jugada con la respuesta del jugador para establecer una nueva situación. 


INSTRUCCIONES 


Al pulsar la tecla RUN se emite un mensaje orientativo de sugerencia 
para quien inicializa la función RND del BASIC. El valor de inicialización 
debe estar comprendido entre 0 y 65535. La introducción de O comenzará 
siempre una secuencia aleatoria diferente; de cualquier otro modo, la se- 
cuencia estará previamente determinada. A continuación, el programa 
emite la llamada indicativa de H (cara) o T (cruz). Finalmente, anuncia 
quién ganó, muestra la puntuación total y continúa con la siguiente jugada. 
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EJEMPLO 


Enunciado del problema 


El jugador (usted) seguirá primero la estrategia de alternar cuando pier- 
da y de jugar lo mismo cuando gane. A continuación, cuando la computa- 
dora adivine su estrategia, invertirá esta última. 


Ejemplo de ejecución 


INICIADOR?1 

(H)CARA O (T)CRUZ?H 

JUEGO 1,TU GANAS, ACIERTOS=1 
(H)CARA O T(CRUZ)?H : 

JUEGO 2,TU GANAS, ACIERTOS=2 
(H)CARA O (T)CRUZ?H 

JUEGO 3,TU PIERDES, ACIERTOS=1 
(H)CARA O (T)CRUZ?T 

JUEGO 4,TU GANAS, ACIERTOS=2 
(H)CARA O (T)CRUZ?T 

JUEGO 5,TU GANAS, ACIERTOS=3 
(H)CARA O (T)CRUZ?T 
JUEGO 6,TU PIERDES, ACIERTOS=2 
(H)CARA O (T)CRUZ?H 

JUEGO 7,TU PIERDES, ACIERTOS=1 
(H)CARA O (T)CRUZ?T 

JUEGO 8,TU PIERDES, ACIERTOS=0 
(H) CARA O (T)CRUZ?H 

JUEGO 9,TU PIERDES, ACIERTOS=-1 
(H)CARA O (T)CRUZ?2T 

JUEGO 10,TU PIERDES, ACIERTOS=-2 
(H)CARA O (T)CRUZ?T 

JUEGO 11,TU GANAS, ACIERTOS=-1 
(H)CARA O (T)CRUZ?H 

JUEGO 12,TU GANAS, ACIERTOS=0 
(H)CARA O (T)CRUZ?T 

JUEGO 13,TU PIERDES, ACIERTOS=-1 
(H)CARA O (T)CRUZ?T 

JUEGO 14,TU PIERDES, ACIERTOS=-2 
(H)CARA O (T)CRUZ?T 

JUEGO 15,TU PIERDES, ACIERTOS=-3 
(H)CARA O (T)CRUZ?T 

JUEGO 16,TU PIERDES, ACIERTOS=-4 
(H)CARA O (T)CRUZ?T 

JUEGO 17,TU PIERDES, ACIERTOS=-5 
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JUEGO 18,TU GANAS, ACIERTOS=-4 
(H)CARA O (T)CRUZ?H 

JUEGO 19,TU PIERDES, ACIERTOS=-5 
(H)CARA O (T)CRUZ?T 

JUEGO 20,TU PIERDES, ACIERTOS=-6 
(H)CARA O (T)CRUZ?H 

JUEGO 21,TU GANAS, ACIERTOS=-5 
(H)CARA O (T)CRUZ?H 

JUEGO 22,TU PIERDES, ACIERTOS=-6 
(H)CARA O (T)CRUZ?T 

JUEGO 23,TU PIERDES, ACIERTOS=-7 
(H)CARA O (T)CRUZ?H 

JUEGO 24,TU PIERDES, ACIERTOS=-8 
(H)CARA O (T)CRUZ? 


Discusión de resultados 


Sigue su primera estrategia hasta el juego 10. La computadora está ga- 
nando cada juego, comenzando con el juego 6. Cuando cambia de estrate- 
gla en el juego 11, ganará usted primero pero la computadora vuelve a 
adivinar su estrategia comenzando con el juego 13. Usted volverá a su anti- 
gua estrategia en el juego 18. Después de unas pocas jugadas, la computa- 
dora vuelve a descubrir su estrategia y después de 24 juegos usted llegará a 
tener un conteo negativo de 8 puntos. 


OBSERVACIONES SOBRE LA PROGRAMACION 


Los resultados de diversos juegos se almacenarán en la variable dimen- 
sionada A. La comprobación de la entrada correcta (H) o (T) se realiza en 
la línea 200. El programa requiere el modo SLOW (lento) de modo que los 
resultados puedan visualizarse continuamente durante la secuencia de IN- 
KEYS. 


LISTADO DEL PROGRAMA 


10: CELS 

15 DIM A(8) 

18 SLOW 

20 PRINT "INICIADOR?"; 
30 INPUT X 

40 PRINT X 

45 RAND (X) 

50 LET J=0 
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LET I=0 

LEFT I=1+1 

IF 1<3 THEN GOTO 140 

LET L=INT ((M*N+1)/2) 

LET W=4*Q+2*L+P+1 

IF A(W)=0 THEN GOTO 140 

LET S=M*A(W) 

GOTO 170 

LET S=2*INT (2*RND)-1 

PRINT "(H)CARA O (T)CRUZ?"; 
IF INKEY$="" THEN GOTO 180 
LET T$=INKEY$ 

IF INKEY$=T$ THEN GOTO 190 
IF T$<>"T" AND T$<>"H" THEN GOTO 180 
PRINT T$ 

LET U=1 

IF T$="T" THEN LET U=-7 

LET J=J-S*U 

LET Z=2-1 

IF. Z=0 THEN GOTO 320 

IF Z>0 THEN GOTO 340 

IF A(W)=0 THEN GOTO 310 

IF S=U THEN GOTO 320 

LET A(W)=0 

GOTO 320 

LET A(W)=M*U 

LET N=M5M 

LET Q=P 

LET M=U 

LET P=INT:*((S*U+1)/2) 

LET A$="GANAS," 

IF:-S=Y4 THEN- LET.'ASf="PTERDES,”" 
PICENT “JUEGO: "SESPSTU- SAS PS] 
GOTO 70 


24. Tiradas 


de dados múltiples 


DESCRIPCION DEL PROGRAMA 


Muchos juegos (v. g. el popular Y AHTZEE, e incluso simulaciones de 
modelos de computadoras están basados en el resultado de una tirada de 
datos múltiples. Este programa permite una tirada de varios dados a la vez 
y visualiza el resultado de cada tirada múltiple: El generador pseudo- 
aleatorio RND del BASIC del Timex/Sinclair 1000 en la línea 90, que 
hace la «tirada», trabaja bien en esta aplicación; sin embargo, si tiene su 
propio favorito, no vacile en sustituir esa línea. 


INSTRUCCIONES 


Al pulsar al tecla RUN se visualiza CUANTOS DADOS para la tirada 
correspondiente. A continuación se introduce un número n (uno o más), el 
programa visualiza n dígitos, en donde cada dígito está comprendido entre 
l y 6 y corresponde a los puntos en un dado. A continuación, aparecerá en 
la pantalla la visualización antes indicada, CUANTOS DADOS, para otra 
tirada. 


EJEMPLO 


Enunciado del problema 


Mantener la tirada de dados, tal como se indica en el listado. 
Ejemplo de ejecución 
DADOS MULTIPLES 


CUANTOS DATOS?5 
hsldals ds 9 

CUANTOS DADOS?3 
6,6,1 

CUANTOS DADOS? 1 

1 

CUANTOS DADOS?1,5 
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CUANTOS DADOS?0 
CUANTOS DADOS?10 
3,4,3,6,H%2.MM.6. 60 
CUANTOS DADOS?2 

6,6 

CUANTOS DATOS?3 
2,5,4 | 
CUANTOS DADOS?7 
3.5461530,2,4 
CUANTOS DADOS?5 
6,3,6,2,5 

CUANTOS DADOS?6 
5172,3,3,3 


Discusión de resultados 


Los resultados serán objeto de visualización. Cuando pedimos 0 ó 1,5 
dados, puesto que ambos números no son admisibles, reaparecerá el men- 
saje CUANTOS DADOS. 


OBSERVACIONES PARA LA PROGRAMACION 


La comprobación de errores de la introducción se realiza en la línea 70. 
El programa es sencillo y corto. La línea 90 indica cómo transformar el nú- 
mero aleatorio generado por RND desde el margen 0-1 a la gama de 
enteros entre 1 y 6. 


LISTADO DEL PROGRAMA 


10. CLS 

20 PRINT "DADOS MULTIPLES" 
30 PRINT 

40 PRINT “CUANTOS DADOS?"; 
50 INPUT A 

60 PRINT A 

70 IF A<>INT A OR A<1 THEN GOTO 40 
80 FOR I=1 TO A 

90 PRINT INT (6*RND)+1; 
100 IF I=A THEN GOTO 120 
TO: PRENT 
120 NEXT 1 
130 GOTO 30 


25. Cuenta corriente 
bancaria de dos usuarios 


DESCRIPCION DEL PROGRAMA 


Se trata de un programa para el control financiero de una familia que 
tiene una cuenta corriente con dos titulares. Actualmente suele darse el 
caso de que trabajen los dos esposos y se plantea, con frecuencia, el proble- 
ma de cómo controlar una cuenta corriente conjunta. Ambas partes pue- 
den disponer de talonarios individuales en los que registran por separado 
todas sus transacciones tales como depósitos, cheques extendidos contra la 
corriente y transferencias entre los dos esposos. Cuando el banco envía un 
estado de cuenta a finales de mes, el siguiente programa ayudará a saldar 
los dos talonarios. Las ecuaciones utilizadas son las siguientes: 


Saldo esposo n.* 2 =saldo mes anterior del esposo n.* 2 + depósitos del es- 
poso n.* 2 - cheques compensados del esposo n.” 2 + transferencias del 
esposo n.* 1 al n.* 2, también 


Saldo esposo n.* 2=saldo bancario corriente-saldo talonario n. 
l —cheques no compensados n.” | — interés devengado (para cuentas 
AHORA) + gastos bancarios. 


El saldo del esposo n.* 2 se calcula por el programa de dos formas y se 
comparan los resultados. Si son los mismos, la cuenta estará saldada. Si son 
diferentes, habrá un error en el estado bancario, o más probablemente, en 
uno de los asientos del talonario. 


INSTRUCCIONES 


Al pulsar al tecla RUN se emite una llamada orientativa correspon- 
diente a las órdenes siguientes: 


Z — Borra todos los registros, comienza programa 

X — Correspondiente al saldo anterior del esposo n.? 2 

X — Correspondiente a cada depósito para la cuenta conjunta por el espo- 
so n.* 2 

V — Correspondiente a cada cheque compensado por el esposo n.* 2 

B — Correspondiente a cada transferencia interna desde la cuenta del espo- 
so l a la del n.* 2; utilizar valores negativos para las transferencias des- 
de el n.*2 al n.* 1 
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N — Correspondiente al saldo del talonario corriente del esposo n.” 1 
M — Correspondiente a cheques no compensados extendidos por el esposo 
azi 

A — Correspondiente al saldo bancario corriente 

S — Correspondiente al interés devengado (cuentas corrientes bancarias de 
AHORA) 

K — Correspondiente a los gastos bancarios, si los hubiera 

D — Visualiza el saldo del esposo n.” 2 calculado en ambas formas. 


EJEMPLO 


Enunciado del problema 


El estado de la cuenta bancaria mensual acaba de llegar y muestra un 
saldo en la cuenta conjunta de $1234.45. El esposo n.* 1 tiene un saldo del 
talonario de $751.23 y dos cheques extendidos por el n.” 1 por las cantida- 
des $15.25 y $10.00 no estaban todavía compensados. El banco abonó en 
la cuenta unos intereses de $6.75. El esposo n.* 2 extendió, durante ese 
mes, tres cheques por las cuantías de $250.00, $125.00 y $50.00, pero el 
último cheque no estaba todavía compensado. El esposo n.” 2 tenía depósi- 
tos abonados en cuenta de $150.00 y $175.00 y una transferencia interna 
del esposo n.” 1 por la cuantía de $300.00. El saldo del talonario anterior 
del esposo n.* 2 era de $203.22 y los gastos bancarios llegaron a ser exacta- 
mente de $2.00. Calcular el saldo del talonario corriente del esposo n.* 2 y 
ver si la cuenta está saldada. 


Ejemplo de ejecución 
CUENTA CORRIENTE DE 2 PERSONAS 


INTRO:+Z ¿XEM BN SM ASS GR 0 DRA 
SALDO ANTERIOR NO0.2?203.22 


INTROD.Z,X,C,V,B,N,M,A,S,K O D?C 
DEPOSITOS NO. 2?150,TOT=150 


INTROD.Z.X+ 0513 BN MA: K, 0, 07: 
DEPOSITOS NO. 2?175,TOT=325 


INTROD.Z,X¡C,V.,B,N,M,A,S,K O D?V 
CHEQUES COMPENSADOS NO. 2?250,T0T=250 


INTRO.Z,X,C,V,B,N,M,A,S,K 0 D2V 
CHEQUES COMPENSADOS NO. 2?125,T0T=375 
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K O D?B 


ENTROD:Z 4 595 
1 00,TOT=300 


>V,B,N,M,A 
TRANSFERENCIAS 1 A 2?3 
INTROD.Z,X 


. )M,A,S,K O D?N 
SALDO NO. 17 


Cs 
27 


ENTROD. ZA 506 


CHEQUES PENDI 15.20 T0T=15228 


A,S,K 0 D?M 
17? 
INTROD,Z.X,C, 


V,B,N,M 
CHEQUES PENDIENTES N 


JA,S,K O D?M 

0. 1?210,TO0T=25.25 
INTROD.Z,X,C,V,B,N,M,A,S,K O D?A 
SALDO BANCARI0O?1234.45 


INTROD.Z,X,C,V,B,N,M,A,S,K O D?S 
INTERES?6.75 


INTROD.Z,X,C,V,B,N,M,¿A,S,K.0 D?K 
CARGAS BANCARIAS?2,TOT=2 


INTROD..Z,Xs EV. BsNM,ASSsK 0: D?D 
SAEDO NO. 2=453.22,493.22 
CUENTA SALDADA 


Discusión de los resultados 


La cuenta está afortunadamente saldada; ambos cálculos llevan al mis- 
mo resultado. Obsérvese que el cheque de $50 extendido por el esposo n.* 
2, que no estaba todavía compensado, nose tuvo que introducir. 


OBSERVACIONES SOBRE LA PROGRAMACION 


La bifurcación se realiza a través de la función INKEY'S y luego, la en- 
trada se compara con las órdenes reconocidas Z, X, C... La línea 580 se 
refiere al error de redondeo del Timex/Sinclair 1000. Obsérvese que en la 
subrutina, que comienza en la línea 1000, todas las variables están puestas 
a 0. ¿Por qué no utilizan la orden CLEAR en su lugar? La razón es que las 
sentencias tales como las que están en la línea 190 o 250 proporcionarían 
el error 2 (variable no definida). 
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LISTADO DEL PROGRAMA 
10-LLS 
15 PRINT "CUENTA CORRIENTE DE 2 PERSONAS” 
20 SLOW 
30 GOSUB 1000 
40 PRINT 


45 PRINT “INTROD.Z,X,C,V,B,N,M,A,S,K,0,D?"; 
50 IF INKEY$="" THEN GOTO 50 

55 PAUSE 5 

60 LET A$=INKEYS 

65 IF A$=INKEY$ THEN GOTO 65 

70 PRINT A$ 

90 IF A$<>"Z" THEN GOTO 110 

95 GOSUB 1000 


100 PRINT "SISTEMA COMPENSADO" 

105 GOTO 40 

110 IF A$<>"X" THEN GOTO 160 

120 PRINT "SALDO ANTERIOR NO. 22"; 
130 INPUT X 

140 PRINT X 

150 GOTO 40 

160 IF A$<>"C" THEN GOTO 220 

UFO. .PRINT "DEPOSITOS: NOS .2t2 os 

180 INPUT P 


190 LET C=C+P 

200 PRINT. PSSS T0T= "56 

210 GOTO 40 

220 IF A$<>"V" THEN GOTO 280 
230 PRINT "CHEQUES COMPENSADOS NO. 2?"; 
240 INPUT Y 

250 LET V=V+Y 

260 PRINT Y;",TOT=";V 

270 GOTO 40 

280 IF A$<>"B" THEN GOTO 340 
290 PRINT "TRANSFIERE 1 A 2?"; 
300 INPUT L 

310 LET B=B+L 

320 PRINT L;",TOT=";B 

330 GOTO 40 

340 IF A$<>"N" THEN GOTO 390 
350. PRINT "SALDO NO..T?*3 

360 INPUT N 

370 PRINT N 

380 GOTO 40 


390 
400 
410 
420 
430 
440 
450 
455 
460 
470 
480 
490 
500 
510 
520 
530 
540 
950 
560 
570 
580 


590 
610 
620 
630 
640 
650 
655 
660 
670 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
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IF A$<>"M" THEN GOTO 450 
PRINT "CHEQUES PENDIENTES N0O.1?"; 
INPUT J 

LET M=M+J 

PRINT J;",TOT=";M 

GOTO 40 

IF A$<>"A" THEN GOTO 490 
PRINT "SALDO BANCARIO?"; 
INPUT A 

PRINT A 

GOTO 40 

IF A$<>"S" THEN GOTO 540 
PRINT "INTERES?"; 

INPUT S 

PRINT S 

GOTO 40 

IF A$<>"D" THEN GOTO 610 

LET. F=X+*C=V+B 

LET G=A-N-M-S+K 

PRINT "SALDO NO. 2="3F;3",";6G 
IF ABS (F-G)<.001 THEN PRINT 
"CUENTA SALDADA" 

GOTO 40 | 

IF A$<>"K" THEN GOTO 660 
PRINT "CARGAS BANCARIAS?"; 
INPUT F 

LET K=K+F 

PRINT ES" TOL="5K 

GOTO 40 

PRINT "REPEAT" 

GOTO 40 

LEE 


Oo0o0o0o0o0o0Oo0Oo0Oo 
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APENDICE 


Técnicas y artificios 
de programación 


Al escribir y depurar los programas de este libro utilicé varias técnicas y 
artificios de programación que desearía compartir con otros usarios del Ti- 
mex/Sinclair 1000 para evitarles algún desconsuelo y frustración. Se citan 
a continuación: 


l. Redondeo 


Si trabaja con números de coma flotante no tiene una decisión sobre 
igualdades de dos números (v. g., IFX = YT HEN...). El TS 1000 trata el re- 
dondeo de forma extraña; v.g. 


100 LET X =3.25 - 3 -.25 
110 PRINT X 


visualizará 1.164... E- 10. Por consiguiente, ha de tomarse siempre el va- 
lor absoluto de la diferencia y compararle con un número pequeño. 


2. Funcionamiento fuera de memoria 


Cuando escriba un programa largo que haga que el Timex/Sinclair 
1000 funcione fuera de memoria, la computadora comenzará a «agonizan» 
y le hará la vida imposible. Uno de. los primeros síntomas de la «escasez» 
de memoria es cuando la computadora ignora la función EDIT (Desplaza- 
do, «SHIFT», 1). Sin embargo, cuando el Timex/Sinclair 1000 comparte la 
visualización y la memoria de programa puede resolver este problema. Por 
ejemplo, si quiere efectuar la edición correctora de la línea 350 mientras 
está trabajando fuera de memoria, haga lo siguiente: 


LIST 350 El listado aparece en la pantalla con la línea 350 en la parte su- 


perior 
CES La pantalla se pone en blanco 
EDIT Inmediatamente después, la línea 350 aparece en la parte infe- 


rior de la pantalla preparada para su edición correctora 
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3. Bifurcaciones múltiples 


El 1000 carece de una GOTO calculada. Para subsanar esta deficiencia, 
numere las diversas bifurcaciones de su pograma en múltiplos de 100 y uti- 
lice la sentencia 


GOTO 100* A 
La variable A debe ser, por supuesto, un entero. 


4. Salida de un lazo 


Si se queda «atascado» durante la ejecución en un lazo «apretado» con 
una sentencia INPUT en el mismo, quizás no pueda salir del lazo sin tener 
que recurrir a desenchufar la computadora. BREAK se interpretará por 
INPUT como un espacio. Si la entrada prevista es numérica, todavía puede 
salirse del lazo tecleando, por ejemplo, Shift Enter (función) LN 0 y obte- 
niendo un error «A». Si la entrada prevista es alfanumérica, su única 
esperanza es probar a llenar la pantalla y obtener un error «5». 


5. Sentencia INPUT 


Preceder a todas las sentencias INPUT con una llamada orientativa. 
Una sentencia INPUT no produce un signo de interrogación, de forma au- 
tomática, como en algunas otras versiones del BASIC. Siga una sentencia 
INPUT con una PRINT para visualizar lo que acabe de introducir. De no 
ser así, nunca lo sabrá. Una secuencia típica podría ser la siguiente: 


10 PRINT «INTRODUCIR X?»; 
20 INPUT X 
30 PRINT X 


6. Utilización de INKEY$ 


INKEYS$ se utiliza como un medio adecuado para la bifurcación del 
programa. Dentro de un lazo con INKEYS$, la computadora debe estar en 
el modo SLOW (o lento) o el programa debe tener sentencias PAUSE; de 
no ser así, la pantalla se mantendrá en blanco sin ningún mensaje o llama- 
da orientativa que sea visible. En este caso, hay una secuencia «a prueba de 
uso indebido» que ha de detectarse; por ejemplo, la clave «X». Trabajará 
en los dos modos SLOW y FAST (lento y rápido): 


100 PRINT «CUAL ES EL SIGUIENTE?»; 
110 PAUSE 10 

120 IF INKEY$ =*“ ” THEN GOTO 110 
130 LET A$ = INKEY$ 

140 PRINT AS 

150 PAUSE 10 

160 IF A$ = INKEY$ THEN GOTO 150 
170 IF A$ < > “X” THEN GO TO 100 
180 REM “X” SE HA DETECTADO 
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Programa n.” 8. Cálculos del desplazamiento relativo para la programa- 
ción en código máquina. 


Introducir los dígitos hexadecimales A a F en mayúsculas o sustituir las 
líneas 280, 290 y 370 y añadir las líneas 281,291 y 371 como sigue: 


280 IF g < =9 THEN LET e$ = CHRS(g + 48) 

281 IF g > 9 THEN LET e$ = CHR + (g + 55) 

290 IF g < 9 THEN LET w$ = CHR$(z — 16 * g + 48) 

291 IF g > 9 THEN LET w$ = CHRS$(z - 16 * g +55) 

370 IF CODE cS(1-1) < =57 THEN LET r = CODE cS(1 -1) - 48 
371 IF CODE cS(1-1) > 57 THEN LET r = CODE cS(1 -1) - 55 
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